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Abstract. In this paper we study reproducing kernel Hilbert spaces of arbitrary smooth- 

rrt ' ncss on the sphere S'' C M''"'"^. The reproducing kernel is given by an integral representa- 

tion using the truncated power function (x • z — t)^ defined on spherical caps centered 

at z of height t. which reduce to an integral over indicator functions of spherical caps as 

CO . studied in [J. Brauchart, J. Dick, arXiv:1101.4448vl [math.NA], to appear in Proc. Amer. 

^^ ' Math. Soc] ioT /3 ~ 1. This is in analogy to the generalization of the reproducing kernel 

to arbitrary smoothness on the unit cube. 
.^^ ' We show that the reproducing kernel is a sum of a Kampe de Feriet function and the 

►^ . Euclidean distance ||x — y|| of the arguments of the kernel raised to the power of 2/3 — 1 

^~i ' if 2/3 — 1 is not an even integer; otherwise the logarithm of the distance ||x — y|| appears, 

i -p^ . For P G N the Kampe de Feriet function reduces to a polynomial, giving a simple closed 

C^ ' form expression for the reproducing kernel. 

Using this space we can generalize Stolarsky's invariance principle to arbitrary smooth- 
ness. Previously, Warnock's formula, which is the analogue to Stolarsky's invariance 
principle for the unit cube [0, 1]'', has been generalized using similar techniques [J. Dick, 
^ Ann. Mat. Pura. AppL, (4) 187 (2008), no. 3, 385-403] 

in 



o 

(N 



1. Introduction and statement of results 



Let S*^ = {x G W^^^ : X- X = 1} be the unit sphere in M'^"'"^ provided with the normahzed 
Lebesgue surface area measure ad, that is Jg^ da^ = 1. The surface area of S*^ is denoted 
by Ud- For future references we define the constant 

><: ..^ _lu;d-i _l T{{d+l)/2) 

^: "- d Ud d ^ r(rf/2) • 

A spherical cap centred at z G S*^ with "height" t G [—1,1] is the set 

C(z;t):={xG§'^:x-z>t}. 
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2 J. S. BRAUCHART AND J. DICK 

Recently, we [5] used the kernel function 

(2) Kci^,y)= f / lc(x,t)(z)lc;(y,t)(z)dad(z)dt, x,yeS^ 

J-i Js-i 

where 1a is the indicator function of the set A C M'^^-'^, to reprove Stolarsky's invariance 
principle (cf. K. Stolarsky [22]) 

2 

dcrd(z) dt 



AT TV 



(3) i=l fe=l a ^-1 ^h „_i 



TV 



|x-y||dcrd(x)do-d(y) 

using reproducing kernel Hilbert space techniques. In this way we provided the theoretical 
underpinning for the particular role this principle plays in the theory of worst-case error 
of numerical integration rules for functions in a Sobolev space of index {d + l)/2 over the 
sphere S'^ provided with the reproducing kernel i^c(x, y) as observed and utilized in [7]. 
This kernel has the very simple closed form (see [5]) 

(4) Kc(x,y) = l-C,||x-y||, x,yG§'^. 

A generalization of the corresponding reproducing kernel Hilbert space to different smooth- 
ness classes has been studied in [TJ Section 6]. This is done via the decay of the Fourier 
coefficients in the series expansion using spherical harmonics of the reproducing kernel. 
However, no extension of Stolarsky's invariance principle for this generalization is known. 
On the other hand, Stolarsky's invariance principle has an analogue for discrepancy 
defined on the high-dimensional unit cube [0, 1]^ with respect to rectangles anchored at 
the origin 0. This analogue is called Warnock's formula [23], see also [TOl Chapter 2]. 
The reproducing kernel for the unit cube is defined analogously to ([2]), where the spherical 
caps are replaced by rectangles anchored at 0. The reproducing kernel Hilbert space 
which corresponds to this kernel consists of functions with square integrable partial mixed 
derivatives of order 1, see [TUl Chapter 2] for details. In [U] this reproducing kernel was 
generalized to include functions of fractional smoothness. This generalization was achieved 
by replacing the indicator function 1[q^z){x) by (x — z)^ for /3 > 1/2. The truncated power 
function (x — z)^ arises from the Taylor expansion of a /3-times continuously differentiable 
function / : [0, 1] — t- M given by 

The Taylor series can be written in a concise form using a reproducing kernel Hilbert space. 
Namely, if we define a reproducing kernel 

13-1 
r=0 



K(x,y) = 5^xV+ / ix-t)i'\y-tt 
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and the corresponding inner product 

then dS]) can be written as 

f{x) = {f,K{;x)). 

In [9] it was shown that this approach yields a generahzation of the geometric discrepancy 
to a fractional geometric discrepancy. 

The aim of this paper is to obtain analogue results for the sphere S"^ as have been shown 
in the cube case. In particular we prove a generalization of Stolarsky's invariance principle 
for function classes of fractional smoothness considered in [3, Section 6] . This case is not 
covered by the results in [22] • 

As in the cube case, the identities 

(6) lc7(z;t)(x) = l[t,l](x-z) = (x-Z-t)'^, X,ZG§^-l<t< 1, 

where x" denotes the truncated power function motivate the following generalization. Let 
(3 > 1/2. For X, y in E>^ we define the function /C^ : S'' x S'^ — )• M by means of 



(7) 



/C;3(x,y):=/ / (x-z-t)^-^(yz-t)^-Marf(z)dt 



It can be checked that /C/3(x, y) < oo for all x, y G S*^ and all (3 > 1/2. Note that for 
(3 < 1/2 this is not true, since /C^(x, y) > /C/3/(x, y) for (3 < (3' and /Ci/2(x, x) = oo. 

The function /C^ is obviously symmetric, that is, /C/3(x, y) = /C^(y,x). Further, let 
arbitrary ai, . . . , ajv G C and Xi, . . . , x^r G S^ be given. Then we have 



N N „-i ^ N N 



J^^aja^/C;3(xj,Xfc) = / / ^^aja^(t-Xj ■z)(J^"^(t-Xfc-z)(J^^Mcrd(z)dt 
j=i k=i -^-^ "^s" j=i fc=i 



1 

■ 1 J^<i 



N 2 



Xj- ■ ZJ'^ 



dad(z)dt > 0. 



^ aj {t 

k=l 

Thus, the function /C^(x, y), is also positive definite. From [2] it follows that /C/3(x, y) is 
a reproducing kernel for each /3 > 1/2 that uniquely defines a reproducing kernel Hilbert 
space ^^/^(/C^, S*^) of functions endowed with a certain inner product which we denote by 

Let 'H/3:='H/3(/C/3,S'^). We note that the reproducing kernel Hilbert space Tip consists of 
all functions /i , /2 : S'^ — t- M which permit the integral representation 

(8) M^)= f [ 9^{z■,t){^■z-t)l^'dad{z)dt, ^ = 1,2, 

J-i J§d 

for some functions gi,g2 '■ ^'^ x [—1,1] — )> M with gi,g2 G L2(S'^ x [— l,l],o"d) of square- 
integrable functions on S'^ x [—1, 1]. See Appendix Rl for a proof of this statement. The 
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function Qi is called the potential function of fi. Thus the inner product for the functions 
/i,/2 e 'His is given by 

(9) {fij2)ic,= f [ gi{z;t)g2{z;t)dad{z)dt. 

For a function / G l-Lp one has 

ll/lk, = {y^^£l^(z;t)|'da,(z)dt| 



1/2 



This follows from the uniqueness of the inner product (-, ■)/c^ induced by the kernel K,p. 

In the following subsections we describe the results in more details. The proofs can be 
found in Section [2] and the Appendix. 

1.1. The smoothness of the Hilbert space l-Lp{K,p,'&^) and the Sobolev space 

ff(S'^). The Hilbert space L2(S'^,(T(i) of square-integrable real-valued functions on S"^ is 
endowed with the inner product and induced norm 






Lacs'*) ■ 



In the following we define the Hilbert space EI*(§'^). Normalized Gegenbauer polynomials 
Pi'^\t) = C(f"^)/^(t)/C[;^~^^/^(l) are orthogonal on the interval [-1, 1] with respect to the 
weight function (1 — ^2^^/2-1 g^j^^^^ thus, related to the unit sphere S'^ as can be seen from 
the Addition theorem for spherical harmonics 

Z{d,n) 

Y^ Y„,fc(x)Y„,fc(y) = Z(rf,n)Pi'^)(x.y), x.yGS"^. 

fc=i 

The real spherical harmonics Y„ ^ form a complete orthonormal basis 

{Y^^l:k = l,...,Z{d,n);n = 0,l,...} 
of the Hilbert space L2(S'^, ad) of square-integrable functions on S'^. Here 

Let s > and let Ag, A'^, . . . be a sequence of positive real numbers such that 

(10) A'„ X n~'^^ as ri — )■ cxD. 

The Sobolev space ]HI*(§'^) over S*^ with index s can be obtained as the completion of the 
family C°°(S'^) of infinitely differentiable functions on S'^ with respect to the norm 



d,s 



■—yU\f)d,s 
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induced by the inner product 

CO Z{d,n) 

(11) a^W=EE(^n)-7SC 

n=0 fc=l 

where the Laplace-Fourier coefficients are given by 

(12) /S = {fMnl)um = I /(x) Yi';i(x)da,(x). 

Obviously any norm equivalent to || ■ H^ ,, will produce the same space of functions. 

For s > d/2 the Sobolev space E[*(§'^) is a reproducing kernel Hilbert space. If we endow 
ff(S'^) with the inner product {■,-)d,s given by ( ITTll . then the corresponding reproducing 
kernel is given by 

oo Z{d,n) oo 

n=0 fc=l n=0 

The reproducing kernel Kd^s satisfies the reproducing kernel properties 

i^rf,,(-,x) e e^(§'^) and (/,Krf,,(-,x)),,, = /(x). 

This is justified by the fact that point-evaluation is bounded in ]H[^(S'^) for s > d/2 by ( TTO|) . 
Imbedding theorems relate Sobolev spaces EI*(S'^) and classical function spaces C'^'(S'^): 
W(S'^) is a subspace of C"=(§'^) iis> k + d/2. 

For further details we refer the interested reader to [7]. Therein, in Section 6, a particular 
instance of a reproducing kernel Hilbert space with inner product and reproducing kernel 
as above has been studied. The reproducing kernel studied there has the property that it 
can easily be evaluated at given function values, since it can be written as a polynomial 
of X ■ y. However, no extension of Stolarsky's invariance principle seems possible for this 
kernel function. 

The following theorem shows that the reproducing kernel Hilbert space "H/j corresponding 
to the reproducing kernel /C^ coincides with EI^(§'^). 

Theorem 1. Let d > 2 and suppose that [5 > 1/2. Then the reproducing kernel /C^ has 
the representation 

oo 

(13) /C;3(x,y) = ^A„Z(rf,n)Pi'^)(x-y), 

n=0 

where 

\n X n~^* as n ^ oo and s = (3 1 — . 

Thus, the reproducing kernel Hilbert space 'H/3(/C/3,S°') uniquely defined by the positive def- 
inite kernel ([7]) coincides with the Sobolev space M.^(E>'^) with smoothness s = f3 — l/2 + d/2 
and the corresponding norms \\ ■ \\d^s O'^d || ■ ||/c„ are equivalent. 

Following [5], we derive a generalization of Stolarsky's invariance principle in the follow- 
ing. 



J. S. BRAUCHART AND J. DICK 



1.2. The Worst-Case Error of Numerical Integration. Let d > 2 and (3 > 1/2. In 
the foUowing we assume that the weights Aq, Ai, . . . are defined by the expansion flT^ . 
The worst-case error of the equal-weight numerical integration rule 



1 ^ 



i=i 



for functions / in the unit ball in the Sobolev space His is defined as 



wce('H^,Q^):=sup 



{|^p'-'-I 



/(x) do-d(x) 



:/G-H^,||/lk, <1 



Invoking the integral representation (^^ of the reproducing kernel /C/j we have 






(x ■ z - t) 



/3-1 



^ £(-.--<-' 



j=i 



do-d(z)dt 



and 



/ /C^(x,y)da,(y)= / / (x-z-t)J-^[/ (y ■ z - t)^"^ darf(y)l da,(z) dt. 

Therefore, the "representer" of the error of numerical integration, 

7^('H/3, Q^;x)/c^ 

1 r Af 

= / /(x-z-t)^-^ 1 ^(x^..z-t)^-i- /■ (y.z-t)J-^da,(y) 
7-1 7§d iV -^ jgd 



da-d(z)dt. 



is of the type of function given in ([8]). The function is called representer since 



(14) 



1 ^ f 



It is well known that for / G 7//3 the following Koksma-Hlawka type inequality holds 



1 ^ r 



[(/, 7^(7^/3, Qtv; ■)/c^)/c^ I < ll/ll/c« [[^('^/S'QAf; Ok: 



pWjc,} 



where equality is assumed for the "representer" of the error of numerical integration of 
the rule Q^ for functions in "H^. See [15] or [lOl Chapter 2] for numerical integration in 
reproducing kernel Hilbert spaces. 
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It follows that 
(15) [vceiHp, Qjv)]^ = ||7^(?//3, Q^; ■)kp\\\^ 

N 



/C/3(x,y) d(Td(x)d(Td(y) - ^y^ ^ib{^,^j) darf(x) 



N N 



j=l k=l 

Since the reproducing kernel /C^(x, y) is a zonal function (that is, only depends on x ■ y), 
see Theorems [H] and [HI the integral /gd/C^(x, y) dcrrf(x) does not depend on y G S'^. Hence 

2 1 ^ ^ f f 

(16) [wce('H/3,Qiv)]^ = THy^y^^/3(xj,Xfc) - / / /C;3(x,y) dad(x) darf(y). 

1.3. Generalizing Stolarsky's Invariance Principle. We define local discrepancy func- 
tion with coefficient /3 of the point set Vn = {xi, . . . , xjv} by 

1 ^ r 



For /3 = 1 we have 



(x-z-t)'[ = lc(z;t)(x) = \ Q 



1 ifxGC(z;t), 
otherwise. 



Thus, the local discrepancy with coefficient /3 = 1 is just the local spherical cap discrepancy. 
Note that the integral appearing in the definition of the local discrepancy above does 
not depend on z and t and can be expressed in terms of a Gauss hypergeometric function. 
(For /3 = 1 it is simply (Jd{C{z^t)) as follows from ([6]).) 

Proposition 2. Let d > 2 and (3 > 0. Then 

Using the local discrepancy function, the representer of numerical integration can be 
written as 

7^('H/3, Q^; x)/c^ = / / (x-z-t)^"^ Ap^,/3(z,t)dcrd(z)dt. 
J-i Js-i 

Hence, using the inner product representation of the worst-case error given in flTSj) and the 
inner product formula in ([9]), we arrive at 

(17) wce(?/^,Q^) = ||7^(7^^,Q;v;x)^J^^ = < / l\Ap^Az,t)fdad{z)dt\ . 
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The right-hand side of fITTj) gives rise to the definition of a generahzed discrepancy which 
for /3 = 1 reduces to the spherical cap h2- discrepancy 



DtiVN):-- 



|PjvnC(x;t)| 



A^ 



^d(C(x;t)) 



1/2 



do-rf(x) dt 



As for the spherical cap L2-discrepancy, one could also study the Lp-norm analogue 



1 JS'* 



|Ap^,^(z,t)|Pda,(z)dt 



i/p 



with the obvious modifications for p = oo. Note that this expression is defined for all 
1 < p < oo if /3 > 1 and all 1 < p < (1 - (5)-^ for < /3 < 1. Let 1 < g < oo with 
l/j> + 1/g = 1. Equation fITK]) is the worst-case integration error for functions with finite 
norm 



K.I3A- 



(?(z,t)|^d(Trf(z)dtl 

1 Js<* J 



l/<? 



with the obvious modifications for q = oo. Here, g is the potential function of /. This 
result can be obtained by applying Holder's inequality to (JH]). 

The combination of the relations f lT6|) and flT7|) gives the generalized invariance principle. 

Theorem 3. Let d > 2 and f3 > 1/2. Then for any N -points xi, . . . , x^v on S"^ one has 



(19) 



J-lJ§d iV J§d 



d(Td(z)dt 



N N 



iv^EE^ 



lX,-,XfeJ 



j=l fc=l 



/C^(x,y) dad(x)d(T<i(y). 



The double integral over the kernel JCfs in ( TT9|) can be reduced as follows. In the following 
let pFq denote the generalized hypergeometric function and (a)^ the Pochhammer symbol. 

Proposition 4. Let d>2 and (3 > 1/2. Then 



/C/3(x,y) dcrd(x)do-d(y) 



-,d-2 



./.-i/2+./2 r((rf + i)/2)r(/3) 



nTiP + d/2) J 7o 



X 



2/3+d-2 



For d = 2 one has 



/C/3(x,y) dcr2(x)d(T2(y) 



A-rf/2,d/2 
^^' /3 + rf/2 '-^ 



22^ 



dx. 



/3 2/3 (2/3 + 1) 
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In general, for d > 2 one has 

/C;3(x,y) do-d(x)dcrrf(y) 



-)d-2 



gd J§d 



2/3 + d-l 



og-i/2+d/2 r((rf + i)/2)r(/3) 

v^r(/3 + ci/2) J 



^ -^1:1:1 



2(3 + d- I 

213 + d 



1 - rf/2, d/2 
f3 + d/2 



1- d/2, d/2 

13 + d/2 



1,1 



For even dimension d the Kampe de Feriet function terminates after finitely many terms. 

If /3 = M , where M is a positive integer, the reproducing kernel /C^ takes on a partic- 
ularly simple form suitable for numerical computations. We recall that this kernel is the 
reproducing kernel for the Sobolev space E[*(§'^) with smoothness s = M + [d — \)/2. 

Proposition 5. Let d >2 and suppose that (3 is a positive integer M . Then 

/CM(x,y) = QM~i(x,y) + {-lY'cd{M) ||x - y||=^"'^^ , x,y e §^ 
where Qjv/-i(x, y) is the following polynomial of degree M — 1 in terms o/x ■ y, 

22M^i (d/2),,,_, A _ M, 1/2 - M, 1 - M l-x^\ 
yM_Ux,yj-2^_^ (rf)2,,_^ '^\3/2-M,2-d/2-2AP 2 J 

and the positive constant Cd{M) is given by 

,_2^,r((c/ + i)/2) r(M)r(M) 



CdiM) 



nTid/2) (1/2),, {d/2),/ 



Moreover, 



QM-i(x,y) dcrrf(x)dcrrf(y) 

. ..M^i JM-iy. fl-M,l/2-M,M + d-l,d/2 
^ ' (3/2),,„/^V l-M,(ci+l)/2,d 



(The 4 F^-function is a polynomial of degree M —1 evaluated at 1.) For d = 2 this relation 
becomes 

/ QA/-i(x,y)da2(x)da2(y)= ,,,,_, ,,2 + ,A.,. ..2 

52 ./s;2 



' ' '' '' M2(2M + 1)' M(2M + 1)M1/2)m ' 
The polynomial part Qm-\ of the kernel Km-, M a positive integer, gets annihilated in 



N N 



Xfe 



/CA/(x,y) do-d(x)dcrd(y) 



i=i fc=i 
if the nodes xi, . . . , x^r form a spherical t-design with t> M — 1. 

A spherical t-design {xi, . . . ,X7v(t)} ^ S'^, introduced in the ground breaking paper [8] 
by Delsarte, Goethals and Seidel, is characterized by 



P(x) darf(x) 



^ Vp(x,^ 



for all polynomials P with degree < t. 
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Seymour and Zaslavsky [21] proved that a spherical t-design always exists if N{t) is suffi- 
ciently large. (From a bound in [S] it follows that N{t) > cf^ for some constant c > 0.) 

Theorem 6 ([6], [IH [121 [131 [H])- Let s > d/2. There exists a constant c{s,d) > 
depending only on s and d such that for any N -point quadrature rule Qn we have 

(20) ^^ < wce(ff , Qn) for allN>l. 
Further, there exists a constant C{s, d) depending only on s and d such that 

(21) wce(ff , Qn) < ^^i^ for all t > 1, 

where the quadrature rule Qn uses an N{t) -point spherical t- design as quadrature points. 

Remark. Korevaar and Meyers [T7] conjectured that to every t there exists a spherical 
t-design with N{t) < Cf^ (C > some constant) points. (A. Bondarenko, D. Radchenko, 
and M. Viazovska propose in [1] a proof of this conjecture.) In light of this the bounds in 
(!20|) and (!2T|) are best possible. 

As a corollary to (ITTl) . Theorem [31 Proposition [SI and Theorem [S] we have the following 
result. 

Corollary 7. Let d>2 and (5 > 1/2. There exists a constant c'(/3,(i) > independent of 
P and d such that for any N -point configuration Vn we have 

c'{P,d) ' '-' <■ ^"' 



Ni/2M^~mM ^ (/^/jAp^,,(z,t)rda,(z)dt) for all N > I. 

Further, there exists a constant C'{(3,d) > depending only on [i and d such that 

where the point set Xn is an N{t) -point spherical t-design. 

Moreover, if 13 = M, M a positive integer, then for any spherical t-design Xn 
{xi, . . . , xat} on S'^ with t > M — 1 and N = N{t) points, one has 



|Ax^,^(z,t)| dad{z)dt 

N N 



-ir'-'ciM) 



|2M-1 



_^ ^^ j=i fc=i 

1.4. Representations of the reproducing kernel. Our ffist result in this subsection 
concerns a closed form representation of the reproducing kernel /C/3 for arbitrary /3 > 1/2. 
Essentially, it can be expressed as a sum of a constant multiple of a (conditionally positive 
definite) (signed) power of the Euclidean distance and a certain integral of a hypergeometric 
function over the sphere S*^ which turns out to be a special case of a Kampe de Feriet 
function. The latter simplifies to a polynomial of the point distance (or their inner product) 
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if /3 is a positive integer, as stated in Proposition [51 The exceptional cases when /3 equals 
L + 1/2 for some positive integer L, can be obtained by a limit process from the general 
result. In this way a logarithmic term in the Euclidean distance is introduced. 

A Kampe de Feriet function in two variables is defined by the double series expansion 



(22) 



pp:r-;/x f '^'1 ! • • • 5 ^-p I Ci , . . . , C^ 

^<i-s;'y[bi,...,bg\ di,...,ds 



111 • • • 1 JlJ. 



00 00 



EE 



(«l) 



m+n 



P'm+n 



x,y 



(Cr)^ Ul)n ■ ■ ■ iU)n ^"^?/" 



where {a)m denotes the Pochhammer symbol. It is a generalization of the Gauss hyper- 
geometric function 



^F<V'^)=E 



The second part of the theorem(s) concerns the ultraspherical expansion in terms of 
normalized Gegenbauer polynomials P„ (t) = Cl^ (t)/ Cl. ~ (1). 

Theorem 8 (general case). Let d > 2. Let 13 > 1/2 with 2/3 — 1 not an even integer > 0; 
that is, (3 — 1/2 = L + e for some integer L >0 and some real < £ < 1. 
The kernel ([7]) has the following representation 



(23) 
where 



Q/3-i(x,y 



/C;3(x,y) = Q/3_i(x,y) + (-1)^+1q(/3) ||x - y 

1 p A/2-/3,l-/3 / 
2/3-1 2-^ M 3/2-/3 'I 



1 2/3-1 



x,yG§^ 



xz— yz 

2— x-z— y-z 



(l-|x-z-|yz) 



1-2/3 



p2:0;l 
■'^^ 1:0:1 



1/2-/3,1-/3 
(t/ + l)/2 



1/2 
3/2-/3 



do-d(z) 



x-y 1-x-y 



x,yG§^ 



and q(/3) zs i/ie positive constant 
q(/3): = 



2-2/3 



r((rf + i)/2)r(/3)r(/3) 



sin(7r(/3 - 1/2 - L)) T{/3 + d/2) r(/3 + 1/2) ' 

For /3 = 1 i/ie kernel ( 123|) reduces to (j4]). 

Lei^ s 6e defined by means of 2f3 — 1 = 2s — d. The ultraspherical expansions 



(24) 



Q^_i(x,y) = 5^AfHs,rf)Z(rf,A;)pf(x-y), 



fc=0 



reduces to a polynomial of degree M — 1 if /3 = M (M a positive integer) and the coefficients 
A\, {s,d) have the following asymptotic behavior otherwise: 



A^\ 



s, d) 



-)d-2 



r((rf + i)/2)r(s + i/2)r(2g) 
r((c/ + i)/2-s)r(s + i) 



r(^/2) 

V{d/2 + 2s; 



k 



-is 



as A; —)■ oo. 
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Thus, the series expansion flM|) converges for all x, y G S*^ if s > d/2. 

Various forms of the coefficients A^ (s, d) can be found in the Proof of Theorem [HI 
The following potential-theoretical constant (also appearing at the right-hand side of 
Stolarsky's invariance principle ([3]) if A = 1) plays a role: 

y§dj§d V7rr(rf + A/2) 

(Its computation can be found, for example, in S. Hubbert and B. J. C. Baxter [16]; also 
cf. Lemma[IIl For the connection to potential-theory we refer the reader to G. Bjorck [3].) 
The expansion of the (signed) power of the Euclidean distance is well-known and follows, 
for example, from S. Hubbert and B. J. C. Baxter [TH]. (For details see, for example, [7].) 
Using the notation of Theorem [S] one has 

(26) (-1)^+^ ||x - y 11--'^ = f^ 1^ Z{d, k)Pi'\^ . y), x, y G §^ 

where (cf. ([25])) 



The coefficient ak{s^ d)/Z{d, k) is alternating in sign ioi k < L and positive ioi k > L + 1. 
The coefficients have the following asymptotic behavior 

^^^-(-1) ^--.(S)^^^/2-.)^ asfc^oo. 

For further references we record the two special cases {L = [/3 — 1/2J) 

ro«^ r r. .^ - i r((rf + i)/2)r(2/3 + d/2- 1) 

^ ^ ^^ ' ^ " 2/3 - 1 r(/3 + {d- l)/2) r(/3 + d/2) ' 

(29) ^^i^.-^)-i~^r'c.m''-'^^,.^.^^^^^^ 



The relation f[28]l is valid for each /3 > 0. The relation f[29]) is valid for 2/3 — 1 not an integer 
>0. 

Theorem 9 (exceptional cases). Let d > 2. Let 2/3 — 1 6e the positive even integer 2L. 
The kernel ([7]) has the following representation: for x, y G S"^ 

(30) /C,+V2(x, y) = Ql-i/2(x, y) + (-i)'^'^2rp^^iy|^ 11^ - y f ' 1^ llx - y II ' 
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where the integral representation of Ql-i/2 is given by 

2L-2n 



QL-i/2(x,y) 



Isr^L 1 (l/2-L)^_„ / x-z-y-z 
2 Z^n=l n (L-n)! I 2-x-z-y-z 



'l-^x-z-iy z) 



-2L 



dcTrffz) 



1 (1/2-L)^ ^oo 1 (l/2)„ / x-z-yz 

2 L! Z^n=l n (L+1)^ I 2-x-z-y-z 



2n+2L 



(l-|x-z-|y z) 



-2L 



■d(Trf(z) 



+ 



(-1)^+^ (1/2)^ 
2 22^ L! 



(x-z-y-z 



,2L 



In 



ll'c-yll 



-2-x-z-y 
anc? i/ie series representation of Ql-i/2 is given by 



-)' + V^(l/2)-^(L + l) 
■ z/ 



do-d(z) 



x,y e 



L-l 



1 (l/2-L)„(l/2)„ 






^^L-n ((d+l)/2)„n! 



1/2 -L 1 + x-y 
1/2 -n ' 2 



1 - 1 (1/2 - L),^ Jl/2)„^^ ^^^ (.L - n, 1/2 - L. 1 + X - y 



L 



1/2- L-n ' 2 

L-fc 



+ - 



l(l/2-L)Jl/2)^^(-l)'= (-L), /l-x-y^^Vl+x-y 



E 



2 ((d+l)/2)^L! ^ A; (1/2 - L), V 2 ; V 2 

2 In 2 + ^(L + 1) + ^1){L + {d + l)/2) 



(-1)^+1 (l/2)Jl/2), /1-x-y^^ 



2 ((d+l)/2)^L! V 2 
V;(l/2)-V'(i. + l/2)), x,yG§^ 



Let s be defined by means of 2(5 — 1 = 2s — d; that is s = L + d/2. The coefficients in 
the ultraspherical expansions 



(31) 



Ql-i/2(x, y) = Y. ^t\L + d/2, d) Z{d, k)Plf\^ ■ y) 



fc=0 



have the following asymptotic behavior: 



r(2) 



A):>{L+d/2,d)r^-2" 



T{{d + l)/2) r(L + {d+ l)/2) r(2L + d) 



T{d/2) 



V{l/2-L)V{L + l + d/2) 
Thus, the series expansion fl3T]) converges for all x, y G S'^ if s > d/2. 
The ultraspherical expansion of the function 



T{d/2 + 2L + d) 



k 



-As 



\L+l 



\2L 



x-y|| ln||x-y| 
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is also well-known; cf., for example, S. Hubbert and B. J. C Baxter [16]. (For details see, 
for example, [7].) Indeed, for a positive integer L one has 

t ^\L+l\\ ||2Li II II V-^0fc(2L,a) i\r>{d)( \ ^ ^d 

(-1) l|x-y|| ln\\^-y\\=2^^——Z{d,k)P'^'{^-y), x,yG§, 

with the following coefficients: 

M^ = i-l)'^W^,^iS') = (-i)--^iM^ [^(L + d/2) + 2 log2 - i;{L + d)] 

and for k < L 

where ak{2L, d)/Z{d, k) is given in ( 1271) . and for A; > L 

hi2L,d) ^ 2L+d-2 r((^ + l)/2) TjL + ci/2) r(L + 1) Tjk - L) 
Z{d,k) y/n T{k + L + d)' 

The coefficients have the following asymptotic behavior 

bk{2L,d) ^ ,jr.^,^, Tiid + l)/2) TjL + d/2) r(L + 1) ^^„,, as A; ^ oo. 
Z((i, /c) y/n 

For further references we record that 

^ . . 1 V 1 (1/2 - LUm^ i-l)'^' (l/2),(l/2)^ 

,,,, '^'^'^ ' ' 2^n-L ((rf + l)/2)„n! 2 ((rf + l)/2),L! 

X (^(L + 1) + 7A(L + (rf + l)/2) - ij{l/2) - ij{L + 1/2)) . 

2. Proofs 
First, we prove the following auxiliary results. 
Lemma 10. Let /3 > 0. Then for — 1 < a < 1 and —l<b<l with a ^ b 

Hf,{a,h):=j\a-t)l^\h-tf-Ut 

= - l + mm{a,6} 1 + af l + fef 2F1 1 /« ; 7— T^T • 

p \ J- + P 1 + maxja, o|/ 

If b = a and — 1 < a < 1, then the integral exists if and only if 13 > 1/2. In this case 

Hp{a, a) = i^j—^ (1 + af^~^ . 
Ifa = b = -1, thenHp{-l,-l) = Q for all P > 0. 
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Proof. Without loss of generality we may assume that — l<a<6<l. The change of 
variable (l + a)n = l+t gives 

Hp{a,b)= f (a - t)l'^ {b - t)^-^ dt 

\a-tf-\h-tf-^dt 

-1 

= (1 + a) / [1 + a - (1 + a) uf~^ [1 + b - {1 + a) uf'^du 
Jo 

= (1 + a)^l + br^ j\l- ur^ (l - i±| n) '" d u. 

The integral represents a Gauss hypergeometric function (cf. [U Eq. 15.6.1]); that is 

Hp{a, 6) = (1 + min{a, b]) (1 + af^' (1 + bf^' ^ ,¥, (\~+Y; ^) • 

If min{a, 6} = 6, then a and b need to be interchanged in the result above. This yields the 
general formula for the integral. 

In the case b = a and — 1 < a < 1 one has (if and only if /3 > 1/2) 

Hp{a,a)= [\a - tff-' dt = {I + af^-' !\l - uf'-' du - ^^ + «^''"' 



j^i JO 2/3-1 

Clearly, i7^(— 1, — 1) = for /3 > 0. This gives the special case. D 

Lemma 11. Let'-/ > 1. Then 

Jsd i§d a/vt r(7 + d) 

Proof. The first part follows from the spherical symmetry of the measure ad and the identity 

||x-yf = 2(l-x-y), x,yG§'^. 

Using the Funk-Hecke formula (cf. [H]), we obtain for the second part 

/ (l+x-z)^da,(z) = ^^^ [\l + ty {1 - t'f'-' dt 
The change of variable 2u = 1 + t yields 

f (l+x.z)^da,(z) = 2^+'^"2+i^^^ f\^+''/^-^{l-uf^-'du. 

Js'i ^d Jo 

The last integral is the Beta function B(7 + d/2, d/2). By (^ and P Eq. 5.12.1] 
;i+x-z)^da,(z)=2" 



^7 . . ..^_ od-mr((rf+l)/2)r(7 + ^/2)r(rf/2) 



y/^ T{d/2) r(7 + d) 

Simplification gives the result. D 
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Lemma 12. Let /3 > 0. Let 1 > a > b > —1. If 2/3 — 1 is not an integer > 0, then 



1 



[l + a 



,/3-l 






'l + a 



.213-1 



iFi 



2/3-1 ^^^V 2-2/3 'l + a 
If 2(5 — 1 = m and m a positive even integer, then 



1 _ /3, 1 _ 2/3. a - 6\ r(l-2/3)r(/3) 



r(l-/3) 



(a-b) 



2/3-1 



/3 



'l + a 



,/3-l 



1 + 6) 



/3 m-1 



r(2/3) 



5^(-l)^(l-/3),(m-l-A:)!(^-^ 



r(l-/3)r(2/3) 

(-i)'-r(/3) 



fc=0 

(a - bf 



a — b 



-irr(/3) ,^ ,,2/3-1 /! + & 



l + a 



^ oo 



Wt /«-'' 



fc=0 



k\ 



l + a 



ijik + (3) - ^P{k + 1) 



In 



a — b\ , \2/3-i 

[a — b) 



r(l - /3) r(2/3) V l + a. 
If P is a positive integer n (that is, 2/3 — 1 zs a positive odd integer), then 



-(l + af (l + 6r2Fi(^ ^^/^^ ;- 



n-l 



[i+ar'ii+brj^ji 



(1-^)* 



a — b 



Proof. First, we observe that by the assumptions on a and b there holds 

and 



o<i±^<i 



l + a 



< < 1. 



l + a 



Thus, the conditions on the argument of the hypergeometric function imposed by the 
variable transformation rules used in the following are satisfied. 

Let 2/3 — 1 be not an integer > 0. Then, by [H Eq. 15.8.4] (and after changing to 
regularized Gauss hypergeometric functions) 

f,{a,by.J-il + ar'{l + br,F,(\-^^'-'-^ 



/3 

r(/3 + i) 

(3 

r(/3 + i) 

/3 



+ /3 '1 + 



\/3-l 



l + a)^"^(l + 6)%FW^^- 



1-/3,1. 1 + fe 



+ /3 '1 + 



:i + a) 



:i+^)' 



vr 



sm 



W2/3-1)) 



r(2/3)r(/3) 



2F1 



1-/3,1 a-^ 
2 - 2/3 ' 1 + a 



r(i-/3)r(i) Vi + a 



a — b 



2/3-1 



2F1 



2/3, /3 Q-^ 
2/9 'l + a 
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In the second hypergeometric function the lower parameter coincides with one of the 
upper ones. Thus, it reduces to (cf. [H Eq. 15.4.6]) 



2F1 



2/3, /3. 



1 p/2/3,/3.a-& 



1 



1 



^/3 



1 



2/3'1 + ay r(2/3)"^'V 2/^ ' i + ay T{2/3) \^ 1 + aJ T{2(3)\l + a 

To the first hypergeometric function we apply the last of the transformations [H Eq. 15.8.1] 

1 + 6 



p ,1-/3,1. a -& 
2^1 1 2 



p , 1 - /3, 1 - 2/3. « - ^ 



2(3^l + aJ~\l + aJ '^ \ 2-2/3 '1 + ^^ 
From the reflection formula for the Gamma function (|ll Eq. 5.5.1]) it follows that 

-r(i-2/3)r(2/3). 



71 



n 



sin (7r(2/3 - 1)) sin (7r(l - 2/3)) 

Taking all these observations into account, we arrive at 



//3(a,6) 



r(/3) (1 + a)^-^ (1 + b)^ r(l - 2/3) r(2/3) 
'1-/3,1-2/3 a- 6\ 1 



1 



r(2/3)r(/3) \l + a 



X 2F1 



2-2/3 '1 



r(l-/3) Vi + a 



2/3-1 



r(2/3) VI + a 



r(l - 2/3) (1 + a) 



2/3-1 



p , 1 - /3, 1 - 2/3. a - 6 

2^ll 2 



r(/3) 



2/3 '1 + aJ r(l-/3)Vl + a 



a — b 



2/3-1 



The first part of the lemma follows after changing back to Gauss hypergeometric function. 
Let 2/3 — 1 = 777, and m an even integer > 0. Then (cf. P, Eq. 15.8.10]) 

-, m — 1 



/,(„, 6) = £(i^ (1 + a^' (1 + i,)^ ^ ^,3+^ 5: (l^'dW™ 



A;)! 



/3 

V i+ay 



E 



r(/3)r(2/3) 



(/3),(2/3), /a -6 



fc=0 



k\ 



1 + a 



r(l - /3) r(l) ^ A;!(A; + m)! VI + a 



fc=0 



ln^_^-,/,(fc + l)+7/;(A; + /3) 
1 + a 



m— 1 



:^^"''"''^^^-^E(-a(™-i-*r^ "-' 



r(2/3) 

-i)"'r(/3) 
r(i-/3) 

a ~ b 



(a-b) 



fc=0 

2/3-1 



1 + 6 
1 + a 



/3 00 

E 

fc=0 



(/3).(2/3), 

A;!(A; + m)! \l + a 



X 



In 



1 + a 



^(A; 



V^(^ + /3) . 



Using identities for the Pochhammer symbol, one sees that 
(2/3), _(m + l), _(m + l), _ (m + 1), 



(fc + m)! {k + m)\ (1) 



m+fc 



(l)„(m + l), m! r(2/3)- 
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This yields 

Ma, b) = ^^^ X. (1 - /5)^(- - 1 - ^)' [-TT-a) 

(-irr(/3) , /l + 6V^(/3), /a-6 



r(l-/3)r(2/3) 



"-)^^-'(|^)E^(fT^)(^(^-«-^(*-: 



izir£(^(„_,)«-L^) fi±iV,F//'.i^^ 



r(l - /3) r(2/3) ' ' V l + aj\l + aj "\l'l + a^ 

The right-most hypergeometric function reduces to ((1 + b)/{l + a))"''. The second part 
of the lemma follows. 

Finally, let /3 = n and n a positive integer. Then fn{a,b) reduces to a polynomial of 
degree n — 1 and [H Eq. 15.8.7] can be applied. (Note that after the transformation the 
lower parameter satisfies 2 — 2n = l — n — M with M > 0; cf. [U Eq. 15.2.5].) 

/,(».. 6) = i (1 + a)"-' (1 + «)•■ .F, ( V+";,^ ^) 

= (1 + a)-' (1 + b)" i , /"'•■-' .F, (ij-M; ff^ 



where it is understood that (0)q = 1. This completes the proof. D 

With this preparations we are ready now to proof Theorem [HI 
Proof of Theorem\^ Let /3 > 1/2. The kernel function (jTj) can be recast as 

^/3(x,y)= / if/3(x-z,y ■z)do-rf(z), 

where 

Hp{a,b):= I {a-tf_^^{b-t)l~^dt. 

Special values /C/3(x, ±x). Let x, y G S'' with x = y. Then by Lemma [TO] 

Hp{^ ■ Z, X ■ Z) = ^n_^ (1 + X ■ Z)^^~^ 

and therefore 

/C/3(x,x)=/ if/3(x-z,x- z)dcrd(z) = — -/ (1 + x- z)^'^~Mcrd(z). 

JSd /p — 1 J§d 
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By Lemma [11] and the duplication formula for the gamma function 

2^^+'^-^r((rf+i)/2)r(2/3 + d/2-i) _ 1 r((d+i)/2)r(2/j + rf/2-i) 

^^^'^'~ 213-1 0Fr(2/3 + d-l) ~ 2(3-lV{l3 + {d-l)/2)V{l3 + d/2)' 

Let X, y G S'^ with x = — y . Then by Lemma [TU] 

/7,(x . z, -X ■ z) = i (1 - |x ■ z|f (1 + |x • z|f -^ 2F1 ^1 - ^'1- i" '"^ ■ ^' 



-V.-l^ .,; v^-ri- ^1; 2.1^ 1 + ^ 'l+|x.z| 

is a zonal function that does not depend on the sign of x ■ z. By the Funk-Hecke formula 
K,(x, -X) = 2i^i j[' (1 - .)^1 + uf-^ .F, (lf/^1; i^) (1 - u'^-' d «. 

The change of variables 

l-u 1 - X , 2 , 

X = ^ — , i-e. u = — — , du = - dx, 

1 + u 1 + X (1 + xy 

yields 



1-1 



-' -^^Jo {i+xr^'~' '^'y 1+/3 '"^r^- 

By da Eq. 2.21.22] and ([1]) 

r r. ^ ,2g+,.ia;,-il 1 r(/3 + l)r(l)r(2/3) A,2/3 + d-l,2/3.1\ 
^''^''' ""^"^ u, /3 22/3+rf-ir(2/3-l)r(/3 + l) ' 'V 2/3 + l,/3 + l ' 2^ 

_ r((rf+l)/2) 2/3-1 A, 2/3 + d- 1,2/3 A 
" v^r(rf/2) /3 ' 'V 2/3 + l,/3 + l '2/ 

For d = 2 the 3 F2-hypergeometric function reduces to a 2 Fi function that simplifies 
further as follows (see [H Eq. 15.4.28]) 

_ r((rf+l)/2) 2/3-1 /i,2/3 1\ _ r(3/2) 2/3-1 ^ r(/3 + 1) 
^^ ' ^ " v^ r(rf/2) /3 '^\/3+l'2y'-0Fr(l) /3 ^^1(1) r(/3 + 1/2) 

r(/3) 



r(/3-i/2) 

Integral representation. Let x, y G S'' with x 7^ y. 

We introduce cylinder coordinates with respect to the Pole p = (x — y)/||x — y||. Then 

1-x-y 12(l-x-y) l||x-yf 1 
||x-y|| 2 ||x-y|| 2 ||x-y|| 2 

y ■ X — 1 

l|x-y|| 
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and one can write 



x=(v/r^^x*,i;), x*e§'^-\ 



z = (Vl-u2z*,n), z* G§"'-\-l < n < 1. 

By construction of these cylinder coordinates there holds that x* = y*. It follows that 

X- z = uv + a/1 — v^ Vl — M^X* ■ z*, 



Jd-1 



y ■ z = — Mf + Vl — ^^ Vl — M^x* ■ z*. 

From the definitions of the these inner products there follows that 

X ■ Z(_u_z*) = y ■ Z(„ z*), for all — 1 < n < 1 and z* G 

Consequently 

X ■ z > y ■ z iff nt; > — ni; iff 2nt; > iff n > 0. 

In other words, one has x • z > y ■ z for every z in the upper half sphere {u > 0) and one 
has X ■ z < y ■ z for every z in the lower half sphere {u < 0). Hence, by Lemma [TO] 

Hf}{x ■ z, y ■ z) = - (1 + min{x ■ z, y ■ z}) (1 + x ■ zf~^ (1 + y ■ z)^~'^ 

„/l-/3, 1 l + min{x-z,y -z} 

X 2^*1 ■ 



1 + P ' 1 + max{x ■ z, y ■ z} 

^/3(x-z,yz) = l(l+x.zf (l + yzf"SFi(^\-/^^;i±^^, n < 0. 

By the Funk-Hecke formula (note that z depends on u and z*) 

/C/3(x,y) = -^ / / i:f/3(x-z,y -z) (l -n^) ^~Mcrrf_i(z*) dn. 

I^d J-l J§d-1 

From the definitions of the inner products x-z and y-z it follows that x-Z(_„z.) = y-Z(uz.). 
Hence 

(33) /C^(x,y) = 2^^^ r / Hpi^-z,yz){l-u'f'-'da,.,{z*)du. 

^d Jo >/§<i-l 

We compute the following integral 

hp{u):= I iJ;3(x-z,y ■z)dcrd_i(z*) 
By our previous considerations (recall, that z = z(m,z*)) 
hp{u) = l^-p (1 + x. zf-^ (1 + y • zf 2Fi(\-//3^; [^^) da,_i(z*), n > 0. 
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-X, X G S'^, the integrand above does not depend on z*; 



We observe, that in the case y 
that is 

(34) h^iu 



uf-' (I 



./3 p /i-/3,l.l-« 



M> 0. 



y=-x /3 ^ 

For the transformation of variable in the hypergeometric functions above we have to dis- 
tinguish between three cases: (a) 2/3 — 1 is not a positive integer (/3 > 1/2 by assumption), 
(b) 2/3 — 1 is an even positive integer, and (c) /3 is a positive integer. Further observe that 



1 +y ■ z 
0< < 1 



for < n < 1. 



1 + x- z 

(The lower bound is assumed when y ■ z = —1; that is, if and only if z = (— Vl — f^x*,f ).) 
Case (a). Let 2/3 — 1 be not a positive integer. First, let y = — x; that is, v = 1. 
Application of Lemma [T2] to (1M|) and the quadratic transformation [H Eq. 15.8.1] gives 



hfi{u) 



(35) 



r(l-2/3)r(/3) 

r(i-/3) 
r(i-2/3)r(/3) 






[l + u 



213-1 



2/3-1 

1 



2F1 



2/3,1-/3 J^ 
2-2/3 ' : 



u 



,1/2-/3,1-/3 o 
l^^H 3/2-/3 '"^ 



r(i-/3) ^-"' ' 2/3 

Hence, using the series expansion of hypergeometric functions 



n >0. 



/C/3(x, -x) = 2 






hp{u) 



(l-u^) 



2\d/2-l 



du 



^,^_i r(l-2/3)r(/3) ^a;,„i 



+ 



r(l-/3) a;, io 

^c.,_,^(l/2-/3)„(l-/3)„ 



2/3-1 (l_^2)<^/2-l 



dw 



2/3 



Wd 



E 

n.=0 



(3/2 - /3)„n! 



u 



2ri 



(1-.=) 



2xd/2-l 



du. 



The substitution u = x^^^ in the integrals yields 



(36) 



u 



2a 



(1 - U^) 



2\d/2~l 



du 



X 



a-1/2 



(1-x) 



d/2-1 



dx 



ir(a + l/2)r(rf/2) 
2 r(a + (t/+l)/2) ■ 



Therefore (using ([T] 
/C/3(x, -x) = 



,,„, r(i - 2/3) r(/3) r((rf + 1)/2) r(/3) r(d/2) 
r(i-/3) 0Fr(d/2) r(/3 + t//2) 



+ 



1 r((rf + l)/2) ^ (l /2-/3),(l-/3)„r(n + l/2)r(rf/2) 

(3/2-/3)„n! T{n + {d + l)/2) 



2/3-1 V?r(t//2) 



E 

n=0 



,2/3- 



.i r((rf+i)/2)r(i-2/3)r(/3)r(/3) 
r(i - /3) r(/3 + d/2) 



+ 



2/3 



^E 

n=0 



(1 /2 - /3) Jl - /3), (1/2), 
(3/2-/3)„n! ((c? + l)/2)/ 
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Comparison with fl39|) below and changing to hypergeometric function yields 



/C«(x, -x) 



92/3-1. (o. , 1 p A/2 - /3, 1 - /3, 1/2. \ 



Taking into account W2\i below, equation fl2^ follows. 
Let y 7^ —X and y 7^ x; that is < w < 1. By Lemma 

r(i-2/3)r(/3) /• ,^^ ^ ^^ ,2/3-1 



^/3^^) = rn m / (x-z-y-z)'' dad-i{z*) 

(37) Ul - PJ J§d-i 

+ X.-. 2/3-1 ^"^^l 2-2/3 ^^ i + x-z Jd-^'^-^(^)' ">0- 
Because of our parameterization, the first integral can be evaluated as follows: 

'r(i-2/3)r(/3)_ f ,^_ _ ^_ ,2/3-1 



ViH / pQ_ ON ■= ^_^^^'^~^'^> dad_i(z*) 

/• / \ 2/3-1 

II l|2/3-l I I ^-y \ A ( * 

= ||x-y||'^ / --z dad_i(z 

y§d-i V x-y J 



l|2/3-l [ I \2/3-l 1 / *N 

y|| ^ / (p-z) ^ d(Td_i(z J 



II l|2^-l 2/3-1 

= ||x — y|| ^ u '^ . 

The second integral is, in fact, well-defined for Re/3 > 1/2 and 2 — 2/3 not an integer < 0. 
Moreover, the hypergeometric function in the integrand has an integral representation for 
1/2 < Re/3 < 1 (cf. P Eq. 15.6.1]). Let 1/2 < Re/3 < 1. Thus, by the Funk-Hecke 
formula and the substitution 2x = 1 — r 

u (\ f (l+x-z)^^~^ A -/3, 1-2/3 x■z-y■z ^^ 
^'''^''^'■= L. 2/3-1 '^\ 2-2/3 ' l + x-z )^-^-^^-^ 



r(i-/3)r(i-/3)7s.-i7o 

r(l-2/3) a;rf_2 f' f' t~^ {1 - ty" {1 - r^f~'^^^-' 



r(l-/3)r(l-/3)a;rf_ii_iio (1 + ni; + v/T^^v^T^^r - 2nt;t)' ^^ 
2/3-ir(l-/3)r(l-/3) a;rf„i V 7 

.1 .1 ^(d-l)/2-l^l-/3-l (1 _ ^)(''-l)/2-l (1 _ ^)l-/3-l 

Jo Jo {1 — wx — zt) 



dtdr 
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where we set 



w = , — , z 



1 + uv + Vl- v^ Vl - u"^ 1 + uv + Vl- v^ VT^V^ 

We note that 



\w\ + \z\ = w + z = -^^=^^—^=^= < 1 itl X ■ z < 1. 

1 + uv + V 1 — v"^ V 1 — M^ 

For d > 2 and 1/2 < Re/3 < 1 the right-most integral above represents an Appell F2- 
function (cf. jTl Eq. 16.15.2]). Moreover, invoking ([1]), we see that 

, . ^ 1 r(2-2/3) . . r(rf/2) / , . ^\ 2/3-1 

^-^^") - 2P-lTil-P)Til-pf V^m-l)/2) {' + ^^ + v/T^ v/T^) 

^^ r((c^ - 1)/2) r(i - /3) r((d - 1)/2) r(i -p) ^ a -2/3,(^-1)72, 1-/3 
r(rf-i)r(2-2/3) n d-i, 2-2/3 

Simphfication gives 

By analytic continuation this identity holds for all /3 > 1/2 except /3 is a positive integer. 
(It is assumed that x • z < 1.) 

Letting a = 2\J\ — v^ Vl — ^^, & = 2tii;, A = 1 + nt; + Vl — f ^ Vl — v? and observing 
that 2A — a — 6 = 2, Proposition [13] yields the following representation in terms of an 
Appell F4 function: 



(38) /,,,(.) = 5^F. {"fiiil^ll, (1 - «') (1 - n^) 



2 2 



Putting everything together, we arrive at 

r(i-2/3)r(/3) „ 



|2/3-l 2/3-1 



r(i-/3) 

1 A/2-/3,l-/3 , 2W. 2N ..2..2 

3T^n d/2,3/2-/3 '(1-^)(1-") 



2/3 
Substitution into fl33l) yields 

^/3 (x, y ) = ^/3,i (x, y ) + /C/3,2 (x, y 



f u 
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where (using ([1]) and functional relations for the gamma function) 

ir ( \ c^^d-\^{"^-'^P>)^{P>) u ||2/3-l /"^ 2/3-1/1 2\rf/2-l^ 

/C/3,i(x,y) =2 — ||x-y||'' / n '^ (1 - ^ ) du 

_ c.,,, r(i-2/3)r(/3) _ ,,_! 1 r^ _ ../2-1 ^^ 

" c., r(i-/3) "'' ^" 2 70 ^ ^ 
_ r((rf + 1)/2) r(i - 2/3) r(/3) r(/3) r(d/2) _ ^/^-i 

V^r(rf/2) r(i-/3) r(/3 + d/2) ""^ ^" 
_ r((ci + i)/2)r(i-2/3) r(/3)r(/3) _ 2/3-1 
v/? r(i-/3) r(/3 + rf/2) ""^ ^" 

(39) -c,(/3)||x-y|| , c.(/3).-2 r(/5 + rf/2) ^ 

and (recall (137|) ) 

^d Jo isd-1 2/3-1 V z- zp 1 + x-zy 



X (l-M^j'^/^-l^^^^^^^,^^^^ 



A linear transformation of variable for hypergeometric functions [H Eq. 15.8.1] yields 

(l + x-z)^^~^ fl- (3,1-2(3 xz-yz ^ 
2/3-1 '^\ 2-2/3 ' i + x-z 

_ (1+yz)'^"' A - /3, 1 - 2/3 y^z-x^\ 
" 2/3-1 '^'V 2-2/3 ' i + y.z / 

Invoking the fact that x ■ Z(_m^z*) = y ■ Z(„^z*) for all —1 < m < 1, we arrive at 

^^'^^"'^^-^IiX.-. 2/3-1 ^^^l 2-2/3 ' l + x-z J 
X (l-n^) ''/'"' darf_i(z*)dt. 
Reversing application of the Funk-Hecke formula yields 

(40) 'C,A-^y) = l 2/3-1 '^\ 2-2/3 ' l + x-z J ^"-^^^^^ 

The quadratic transformation of variable [H Eq. 15.8.1] and a change of variable z i— t- — z 
give the symmetric form 

pyi/2-/3,l-/3V x-z-yz ^2 



2l_2/3 r 2-iy 3/2-/3 ' V2-x-z-yz 

2/3 — 1 j§d (2 — X ■ z — y ■ z 



(41) /C;3,2(x,y) = 7^^5 ^/ -^32^ ^d(Td(z). 
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The constant in front of the Euchdean distance in fl39l) can be rewritten as follows. Let 
L = [/3 - 1/2J be the largest integer < f3 - 1/2. Then f3 -1/2 = L + e ioi some < e < 1. 
Application of the reflection formula for the Gamma function gives 

j.(. ,._a. r(/3 + l/2)r(l/2-/3) ^ _^ n ^ (-1)^+^ TT 

^' ^' r(/3 + l/2) r(/3 + 1/2) sin(7r(/3 + 1/2)) r(/3 + 1/2) sin(7r£) 

and therefore 

(-1)^+1 21-2/3 r((d + 1)/2) r(/3) r(/3) 



(42) c,(/3) 



sin(7r(/3 - 1/2 - L)) 2/3 - 1 r(/3 + rf/2) r(/3 - 1/2) ' 



Using the series representation of the Appell F4 function (cf. Proof of Proposition [T^ 
in (138|) . we get 



zp — i Wrf Jo 



2 CUd-l y-^ y-^ (1/2 P)m+7i(^ P)m+n (^ 2\™ 2nt 



2/5-1 ^-^ ;^o^ (rf/2)j3/2-/3)„m!n! 
where 

6^,n:= r«^"(l-«^)"^'/'^'d«. 

The change of variable x = v?' gives a beta integral. Hence 

6„,,„ = - /" x"+i/2-i(l-x)'"+^/'-Mx = -B(n + l/2,m + rf/2) 
_ ir(n + l/2)r(m + rf/2) _ ir(l/2)r(rf/2) (l/2)„(rf/2)„ 



2 r(m + n + (rf + l)/2) 2T{{d+ I) 12) {{d + l)/2)^^„ " 
Taking into account ([I]), we arrive at 

(43) K..(x. y) ^ ^^ E g („^i)/2)_,3/2_^)^„,„, (l - « ) - 

The double series expansion has the structure given in (!22|) . (This gives the first part of 
the theorem (for the case 2/3 — 1 is not an integer > 0).) 
Reordering the sum with respect to constant m + n gives 

lAA^ r I^.A 1 ^^^ (1/2 - /3)„(1 - /i),. [^ f^n^, (1/2),. ^.,„ ,, ,,„-„ 

(44) K«,(x. y) = ^a^T Ij ((d + i)/2)„„! [Ij (,J (37237J);" (1 - " ) 

Using [H Eq.s 15.8.1 and 15.8.7], the square-bracketed equals 

h ^.2xn p /^-n,l/2. t^' \ _ /-n, 1 - /3 2A _(V^),^ p /-ri, 1 - /3 

(1-M ^Fi ^3/2 - /3' ^;^3T J - '^A 3/2-/3 '^ j " (3/2 - /3)/^^ 1/2 - n '^ 



t;2 
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By definition oi v, we obtain 

U^) K fxvl- ^ f^ (1/2 - /3) Jl - a /-n,l-/3. 1-x.y \ 

(45) /C,,2(x, y) - ^^-^ ^ ((rf+l)/2)„n! '^' ^ 3/2 - /3 ' ^2^ J " 

We conclude Case (a) by deriving the expansion of tlie kernel /C/3(x, y) in terms of 
spherical harmonics. Evidently, this kernel is a zonal function as it depends on the inner 
product X ■ y. Thus, the expansion takes on the form of a series in terms of ultraspherical 
polynomials. Since the expansion of the (signed) power of the Euclidean distance is well- 
known, we focus here on the expansion of /C/3,2(x, y). 

By Proposition [16] 

' " (1/2 - /3) Jl - /3)„ afc,„ 

2/3-lf^f:^ ((rf + l)/2)„n! Z{d,k) 



r hr v^ - ^ Y^Y^ (V2-/J)„(l-/J)„ ak,n :7(^ JApid)(,^ ,A 

oo 

= ^Af(.,rf)Z(rf,A:)pf(x-y), 

fc=0 

where (using /3 where it is convenient) 

.(2). .._ 1 v (i/2-/3)Ji-/3). g^ 

' ^ ' ^' 2/3-1 £-^ ((rf + l)/2)„n! Z(d,A:) 



and 



n!(l-/3),(d/2), /A;-n,fc + l-/3,fc + d/2 A 

2-a(^Ur^-fc)!' ^ A: + 3/2-/3,2fc + d 'V' 



Z(rf,fc) (3/2 

Let /cq be the smallest integer with /cq + 1 — /? > 0. Then we can use the integral 
representation of afc„/Z((i, k) given in Proposition [TBI if /c > k^. Let /c > k^. Then 



4(2)r, ,v_ 1 f. (1/2 - /3) Jl - /3), g,,, 
' ^' ^' 2/3 -l£^ ((c/ + l)/2)„n! Z(rf, fc) 



1 r(3/2-/3)r(rf) - (i/2-/3)Jl-/3)^ n\ 



2^-ir(l-/3)r(l/2)r(rf/2)r(d/2)2 ((ci + l)/2)„n! (rf/2),(n-fc)! 



1 /•! 



X / / t'=+l-^-V+"'/2-l(^_^)l/2-l(^_^)fc+d/2-l(^_^^^n-fc^^^^ 



^0 

oo 



1 r(3/2-/3)r(d) 1 ^ (i/2-/3)„(l-/3)„ 

2/3 - 1 r(l - /3) r(l/2) r(rf/2) r(d/2) (d/2), £^ ((d + l)/2)„ (n - A;)! 



Jo 



X / / t^+'-^~\^+'''^-\i-tf^-'{i-xf^"'^\i-xtr-Uxdt. 
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We compute 

f. (1/2 - /J) Jl - /3)„ _ _, ^ ^ (l/2-/3W(l-/3W .. _ ..n 
2^^((d+l)/2Un-ky.^ ""> ^ ((d+l)/2)„^,n! ^ ^^ 

_ (1/2 - /3),{1 - /3), ^ (fc + 1/2 - /3) Jfc + 1 - /3)„ . 

iid+l)/2), £^ (A: + (rf+l)/2)„n! ^ ^ 



(l/2-/3),(l-/3), 



((d+l)/2) ^ 



p /A: + l/2-/3,A: + l-/3 _ \ 
^\ k + {d + l)/2 '^ ^^J 



r7;rnn7^v ^"'^^ '^n k+(d+i)/2 '^ "^*'- 



((rf+l)/2), ^^^^ ''\ k+{d+l)/2 

Hence, for A; + 1 - d/2 - 2/3 > (and k > (3 - 1), we have 
(46) 
A^^\s,d) 



(2) , ^. _ 1 Tid/2 - s) T{d) {d/2 - s),i{d + l)/2 - s). 



2 r((d + l)/2 - s) r(l/2) Tid/2) Tid/2) ((rf + l)/2),(rf/2), 

'^(^ + (^ + l)/2).n!yo io (1-xtr ^- 

Invoking integral representations for the hypergeometric 2 Fi and 3 F2 functions and a 
transformation for 3 F2 functions of unity argument (cf. [191 Eq. 7.4.4.83]), we arrive at 

CI tX/ (J. 6 



JO (1 - xt) 

r(2s)r(fc + ci/2) r,.+i/2-w. ,^1/2-iFr -^'2^ -Adf 

- r(A: + rf/2 + 2s) io ^^^ 2FiU + rf/2 + 2s'^jd^ 

_ r(s + l/2)r(l/2)r(2s)r(A: + rf/2) / -n,s + l/2,2s 

r(s + 1) r(A; + d/2 + 2s) ^ 2 Vs + 1, /c + rf/2 + 2s' 

r(s+l/2)r(l/2)r(2s)r(fc + rf/2) ik + id + l))^ f -n,l-s,l/2 

r(s + 1) r(A; + d/2 + 2s) ik + d/2 + 2s)„ ' ^ 1^1 + s, fc + (d + l)/2' 

Substitution into (146!) yields 
(47) 
Af\s,d) 



(2), ^,_ 1 r(rf)r(d/2-s)r(s + i/2)r(2s) (rf/2-s),((rf + i)/2-s). 



2 r(rf/2) r((rf + l)/2 - s) r(s + l) r(rf/2 + 2s) Hd + l)/2)^id/2 + 2s), 



X 



oo 



n=0 



(^ + 1/2) Js)„ ^ / -n,l-s,l/2 



(A; + d/2 + 2s)„n!' \l + s,k + id + l)/2 



For asymptotic analysis of A^. (s, rf) for large k we reverse to (146|) . The infinite series is 
already in a hybrid-asymptotic form (as k becomes large). For the curly braced expression 
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we can, thus, write 

{s + 1/2) s 



I ■ ■ ■ I = r r f+'/'~'x'^-' (1 - tY^'-' (1 - x)''-'^'-'dxdt + 



k+{d+ l)/2 
T{s + 1/2) r(l/2) r(2s) T{k + d/2) {s + 1/2) s 



n{k) 



r(s + l) r{k + d/2 + 2s) k + {d + l)/2 

where 



n{k), 



^W-2.(i + fc + (rf+l)/2)J, + l),X X (1-xt)— ^ ^'^^'- 



n=0 

Since (1 — xt)""*"^ < 1 and (n + 1)! > n\, we have the estimate 



< n{k) < 



T{s + 1/2) r(l/2) r(2g) r(A; + d/2) 
r(s + l) T{k + d/2 + 2s) 



fl + s + l/2,l + s \ 
^^\l + k + {d+l)/T^J' 



Observe that this hypergeometric function at unity argument decreases as k gets larger. 
For fc > 2s + 1 — d/2 we can apply [U Eq. 15.4.20] to evaluate the hypergeometric function 
at argument unity: 

ci(k) < r(g + 1/2) r(i/2) r(2s) rjk + d/2) T{i + k + {d + 1)/2) r(fc + d/2-i- 2s) 

^ ^>- T{s + 1) T{k + d/2 + 2s) T{k + d/2-s)T{k + {d+l)/2-s) ' 

We conclude that 

^r(. + i/2)r(i/2)r(2.)r(<. + d/2) ^,,,^^ 

r(s + l) T{k + d/2 + 2s) ^ ^ '' 

Hence, 
{d/2 - s),{{d + l)/2 - s), I I r((d + l)/2) Tid/2) Vis + 1/2) r(l/2) r(2s) 



{{d+l)/2),id/2), I f T{d/s-s)Tiid+l)/2-s)T{s + l) 

T{k + d/s - s) T{k + {d+ l)/2 - s) 



T{k + {d+ l)/2) r(A; + d/2 + 2s) 
Putting everything together, we obtain 



{i + o{k-')) 



(2), ^,_ 1 r((d + i)/2)r(rf/2)r(d)r(s + i/2)r(2s)r(s + i/2)r(i/2)r(2s) 



A','{s.d) 



2 T{d/2) T{{d + l)/2 - s) T{s + 1) T{d/2 + 2s) T{{d + l)/2 - s) T{s + 1) 
T{k + d/s - s) T{k + {d + l)/2 - s) 



X 



r(A; + {d + l)/2) r(A; + d/2 + 2s) 
Application of [H Eq. 5.11.12] yields that as /c — )> oo 



{l + 0{k-^)), ask^oo. 



,(2). . 1 r(i/2) T{{d + 1)/2) r(rf) r(s + 1/2) r(s + 1/2) r(2s) r(2s) ^_,, 



Al''{s,d) 



2 r{{d + 1)/2 - s) r((d + 1)/2 - s) r{d/2 + 2s) r(s + 1) r(s + 1) 
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The duplication formula for the gamma function gives 



(4 



.(2) 



Al'\s,d) 



ld-2 



r((rf + i)/2)r(s + i/2)r(2s) 



1 2 



Tid/2) 



T{d/2 + 2s] 



k 



-4s 



T{{d+l)/2-s)T{s + l) 

This completes the proof of Case (a) of the Theorem |8l 

Case (b). This is the content of Theorem HI 

Case (c). Inspection of Case (a) of this proof shows that it is save to take the limit as 
/3 goes to a positive integer M to obtain (cf., in particular, ( 139|) and (142|) ) 



/CM,i(x,y) = lim^/C^,i(x,y) = q(M) ||x - y 



|2M-1 



where (note that L = M — 1] 



\M 



)1-2A/ 



r((rf + i)/2)r(M)r(M) 



sin(7r(M - 1/2 - L)) 2M - 1 T{M + d/2) T{M - 1/2) 
^^T{{d + l)/2)T{M)T{M) 



-1 



r(M + rf/2)r(M-i/2)' 

Similarly (cf. (jUD, (03]) and (gS])) 



^A/,2(x,y) = lmi^/C^,2(x,y) 



)1-2A/ 



A/2-Af,l-M / x.z~y.z 
^^\^ 3/2 -M ' 1,2-x.z-y.z 



2M - 1 jgd 

^ oo oo 

2M - 1 ^ ^ 



vl-2A4' 



dfTd(z) 



m=0 n=0 



(2-x-z-y -ZJ 
(1/2 - M)^^J1 - M)^^„(l/2)„ 
((^ + l)/2Wn(3/2-M)„m!n! ^ " > 



,2n 



(49) 



2M-1 ^^ 

n=0 



(l/2-M)Jl-M)„ /-n, 1 - M. 1-x-y 



((rf+l)/2)„n! 



Fi 



\^ 3/2 - M 



(The last series terminates, since (1 — M)^ = for n > M.) 

Since /Cm,2(x, y) is a polynomial of degree M — 1 and Proposition 
positive integer, we can simply write (cf. Proof of Theorem |8]) 

M-l 

/CM,2(x,y) = J2^k\^,d)Z{d,k)Plf\^-y), 



applies if /3 is a 



fc=0 



where (using /3 where it is convenient) 



M-l 



(50) 
and 



I (2) 



A):>{s,dy. 



2M 

n\il-M),id/2\ 

Z{d,k) {3/2-M),{d),,{n-ky. 



1__ '^^ (l/2-M)Jl-M)„ g,,, 
-1^^ ((rf+l)/2)„n! Z(d,A;) 



^k.n 



p fk-n,k + l- M,k + d/2 \ 
^ \ k + 3/2-M,2k + d 'y- 



30 J. S. BRAUCHART AND J. DICK 

(Note that the 3 F2 hypergeometric function above reduces to a polynomial of degree M — 
l-k.) 

Furthermore, the relations for /C/3(x, x) and /C^(x, — x) also extend to the case when /3 
is a positive integer M. This completes the proof. D 

For the proof of Theorem [9] we need the following auxiliary result. 

Lemma 13. Let d > 2. Let L be a positive integer, n an integer > and e a real number 
with l^l < 1 and e ^ 0. Then for x, y G S*^ with y 7^ x and y 7^ — x there holds 

I xz-yz \ 

^ or o dfftj Z) 

■^ (1 - ix ■ z - iy ■ z) 

_ (V2)n „2n -p /^n - L - £, n - L + 1/2 - £ oA 

-((rf + l)/2)/ '^\ n + {d + l)/2 '^ ^'1 

Proof. We use the same notation as in the Proof of Theorem [HI 

Let (7„(x, y) denote the integral we want to compute. Then applying the Funk-Hecke 
formula twice, we obtain 

(^x-z-^yz ) 
sd (1 - ix ■ z - |y ■ z) 
Ud-i [^ f (vuf"^ (1 - u 



^n(x,y)= / ^ ^ ,2n-2L-2e darf(z) 



/ T. n 9 r. 9 ,2n-2L-2. d^d-l(z )du 

Jsd-i h _ ./I -7,2, /I -7/2x* . Z* 1 



c^d J-i y§d-i ^1 _ v^l -i;Vl -n^x* ■ z*) " 
a;,_i a;,_2 r /" (^^^)'" (1 - u^f'~' (1 - r' 



dr dw. 



The standard substitution 2x = 1 + r gives 

O fx Vl - ^^^O-^-^^^^^ 2n /" ^'" (1 - ^'f'~' 
yny^,y)- '^ ^ / ^ ^ ^ ^.2n-2L-2e 

/•I ^(d-l)/2-l (^ _ ^N(rf-l)/2-l 

X / K — „r „ dx du. 

/ ; s- ; PT \ 2n—2L—2e 

H _ ^v^HZv^EmL-j;^ 



l+v'l--«^Vl-«2 

The inner integral represents (cf. [H Eq. 15.6.1]) 



r((d-l)/2)r((rf-l)/2) /2n-2L-25,(rf-l)/2.^VrHZvlE^A 
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which can be written as follows when using the quadratic transformation [1] Eq. 15.8.13] 

X. T? fn - L ~ e,n - L + 1/2 - s /. 2\ h 2^^ 
Taking into account ([I]), we arrive at 

X 2~ 



T^ fn- L - e,n - L + 1/2 - e /, 2\ /i 2^^ i 
Fi(^ \^/2 ;(l-t; ) (1-M )J dw. 



The change of variable x = 1 — u yields 



„ /ra - L -£,r2 - L + 1/2 -£ / 2\ \ i 

X 2F1I ' ^ -2 ; (1 -17 jx 1 dw. 

This is the integral representation of a 3 F2-hypergeonietric function (cf. [H Eq. 16.5.2]): 

. ._ r((rf+l)/2) ^,„ r(c;/2)r(n+l/2) /n-L-e,n-L + l/2-£,rf/2 . _ A 
^"^ ' ^^ ~ y? r(rf/2) r(n + (d + l)/2) ' ' V d/2, n+{d+ l)/2 ' ' 'J 

Because of an upper parameter coinciding with a lower parameter, the 3 F2-function reduces 
to a 2 Fi-function 

'"^"'^^-((rf + l)/2)/ ^^\ n+(rf+l)/2 'l-'^l 

This completes the proof. D 

Proof of Theorem\^ We use the same notation as in Theorem [8] and its proof. The result is 
obtained by a limit process as 2/3 — 1 tends to an even integer 2L. That is, let /3 — 1/2 = L+e 
with e sufficiently small and e 7^ 0. We consider limits as e — ?■ 0. 
Let y = X. Then (1281) implies that 

r ( \ V r ( \ ^ V{{d+l)/2)V{2L + d/2) 

^L+i/2(x,x) = lim/CL+i/2+e(x,x) - 



,_^o -^^/-^^^ ' ' 2LV{L + d/2)V{L + {d+l)/2) 
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Let y = —X, X. Then fl29l) implies that 



"^ ^ '^^^ 2£^(l/2-/3)„^,((d + l)/2)„n! 

- (--[\L+1 jr , 1 /2 I ^^02i+2e _ ^ V ^ (V2-L-£)^(l/2)^^ 



2^^-^-^ ((rf + l)/2)„n! 



Hence 

^L+i/2(x, -x) = hni;C/3(x, -x) 



e-^O 



-lim [11 (V2-L-£)^(l/2)^ f_ni+VYr I 1 /2 + rl22^+2e 
-l^^o\2e ((ci + l)/2)^L! +^ ^^ c,[L + l/2 + e}2 

If- 1 (l/2-LUl/2)„ 
2t^^-i> ((rf + l)/2)„n! 

The hmit above is the same as in ( 15^ below but for f = 1 (that is, ||x — y|| = 2). We can 
use fl55l) below with v = 1. 



r rx x^- ^V ^ (i/2-L)ji/2)„ 



1 (l/2-L)Jl/2)„ , (-1)^+1 (l/2),(l/2) 



L 



2^^n-L i{d+l)/2)y 2 ((rf+l)/2)^L! 

X ( V(l/2) + ^(i^ + 1/2) - ^(L + 1) - ^(L + {d+ l)/2) 

This shows ( l32l) . 

For the remaining proof we assume the following: let y 7^ — x and y 7^ x; that is 
< y < 1. 

The integrand of the integral representation of Q/3_i in Theorem |H] splits into two non- 
critical and one critical part as follows. 

L_ F A/2 -13,1-13. ,\_ 1 f^ (1/2 - /3) Jl - /3) ,. 
H 3/2-/3 '"^ i 2/3-1^ (3/2-/3) n! 



2/3-1^ \ 3/2-/3 '-; 2/3-1^ (3/2 - /3)„n! 

If- (l-/3)n 2n If- 

2^ (n + 1/2-/3) n! 2^ 

(51) 



2^(n + l/2-/3)n! 2^ („ - L - e) n! 



lv (l/2-g-^). 2n l (l/2-e-L) ,^ 1 f. (1/2 -5 -L) ,, 
2^(L-n + 5)n! 2 (-5) L! 2 ^^^^ (n - L - 5) n! " 
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Consequently, 

1 



H\W^' \ 3/2-/3 '^ 



l {l/2-e-L) ,, 
"2 {-e)L\ 



L-l 



-X 



2^^ (L-n)nl 2 ^^ (n-L)n\ 

n=0 ^ ' n=L+l ^ ' 



■X 



ly4(l/2-L)„^,„ 1 ^ (l/2-L)„^,„ 



E 



E 



1 ^ 1 (1/2 - L)^_„ ^^2^^_. _ 1^,^ |. 1 (1/2 - L)^^^ 



1 ^— \ 1 I 

2^ra" (L-n)! ^"^ ^ 2"' ^^n {n + L)\ 

n=l ^ ' 71=1 ^ ' 



X 



2n 



1 V^ 1 (l/2-L)j^^„ ('r^'l^^" - 1 (V^-L)^ ^^ ^ ^ l^/^A 



-^ n {L — n)\ 



L^2L V^ 1 (V2)„ _^2n_ 



t^(^ + l)n^ 



The critical contribution to the desired limit (as e — > 0) for ICp are due to the singled 
out term above and the second part in (|23l) . 

Integral representation. Hence, by Theorem [H] and (13 9 p 



(52) 
^/3(x,y) 



1 p / 1/2 - /3, 1 - /3. / xz-yz y \ 1 (l^/3)z. / xz yz 

2/3-1 2J^1\ 3/2-/3 ' V2-x-z-yzy ^ "^ 2 (-e) r(L+l) I, 2-x.z-yz 

1 \ 1-2/3 



2L 



(1-ix-z-iyz) 



do-d(z) 



r(i-2/3)r(/3) 21 

+ L1 r(i-/3) [(-^-^-)] 



1 (l-/3)i / X 



/3-1/2 2(-£)r(L+l) l2-x-z-y-z 



>z-y-z 



2L 



(1-i 



1 \ 1-2/3 

x-z- ^y -zj 



d(Td(z) 



The integrand in curly braces above can be rewritten as follows 



2^-^^ (1 - /3)^ (: 



x-z-y -zj 



s2L 



(53) 



Since 



L'- (2 — X ■ z — y ■ z 



2L+1-2/3 



22/3 



L!(-£)r(l-2/3)r(/3) / (x-z-yz 



2 X /3-1/2-i 



r(l-/3)(l-/3)^ I (2-x.z-y.z) 



(l-/3)i 



r(/3) 



(/3). 



r(/3-L) 



-^) r(-L -e) = i-s) r(-5)(-e)„^ = r(l -e)-^ 



-ly 



+ e) 



T{l-e)T{l + e) 
T{L + l + e) 
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, L!(-e)r(l-2/3)r(/3) 

r(i-/3)(i-/3)^ 



L!(-£)r(i/2-/3)r(i-/3)r(/3) 



■1) 



r(i/2)r(i-/3)(i-/3)^ 

^ L\ {-e) r(i/2 - /3) r(i - ^) r(/3) r(^ - l) 

r(i/2)r(i-/3)r(/3) 
^r(L + i)(-5)r(i/2-/3)r(/3-L) 



^ ^ r(i/2) 

r(L + 1) r(i - e) r(i + e) r(l/2 + e) 

r(i/2)r(L + i + £) 
r(L + i) r(i/2 + £)r(i + £)r(i-£) 
T{L + i + e) r(i/2) r(i) r(i) 



Let 



G(x;e): = 



r(L + i) r(i/2 + £)r(i + £)r(i-£) 
Y(L + i + £) r(i/2) r(i) r(i) 



X 



-1/2 <e < 1/2, x> 0. 



For X > this function assumes the value 1 at e = and is differentiable at e = with 



dG{x;e) 



de 



e=0 



-i){L + 1) + ^(1/2) + V'(l) - ^(1) + Inx = Inx + V(l/2) - V^(Iv + 1), 



where 'il){z):=T'{z)/ T{z) denotes the Digamma function. 
Hence, taking the hmit as e — )■ in ( 153|) . we obtain 



hm 

£-S>0 



2 " {l/2 — L)^^_ _ _n2Li, ^^V^-x-z-y-z 



2 L\ 

-1)^+1 (1/2) 



2 22^ L! 
21n ||x — y|| + In ( 



(x ■ z — y ■ z) hm 

£-5>0 



'x-z-y -z) 



G{ 



xz— yz 



x-y 

l|x-y|| 



2 — X ■ z — y ■ z 



+ V'(l/2)-V'(i^ + l) 



provided x ■ z — y ■ z 7^ 0. 
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It follows that after interchanging limit and integration in ( 152|) (cf. Appendix [C]) , we get 

^L+i/2(x,y) = lim/C^(x,y) 



1 y^i 1 (l/2--L)^_„ / x-z-y-z 

2 Z^n=l n (L-n)! \^2-xz-yz 

^2L 



2L-2n 



S'^ (l-|x-z-iy-z) 

1 {1/2-L)l Y-oo 1 (l/2)„ /^ x-z-yz 



2 

Eoo 
n=l n (L+l)„ \2-x-z-y-z 

(l-|x-z-iy- z)"^^ 



dcrd(z) 

2n+2L 



■dad(z) 



\L+1 



:i/2), 



2 22iL! ./§d 
x-y 

In ' 



[x-z-y -zj 



n2L 



l|x-y|| 



2 — X ■ z — y ■ z 



-)' + V^(l/2)-^(L + l) 
■ z/ 



do-d(z) 



+ (-1) 



l)-^(^ln 



22iL! 



2L 



x-y|| / (x-z-y-z)' do-rf(z). 



The right-most integral evaluates as (cf. (EH])) 



2L 



(x-z-y-z) dcrd(z 



(1/2), 



|x-y| 



|2L 



((rf+l)/2)^ 

This completes the proof of the integral representation. 

Series representation. From (jHj) and the following relations we obtain the following 
representation 



Q/3-i(x,y) 



2/3 



n=0 



(l/2-/3)Jl-/3)„ (l/2)„ /_n,l-/3 2 

((rf + l)/2)„n! (3/2-/3)/ ^ U/2 - n ' 



which can be simplified to (using /3 = L + 1/2 + £:) 



^ oo 

Q^„i(x,y) = --5^ 



1 (l/2-L-e)Jl/2) 



n=0 



^_L-e ((d+l)/2)„n! 



/-n,l/2-L-£ 2^ 

'^\ 1/2-n '1-^J- 

This sum can be also split into two non-critical and one critical part regarding a limit 
process /3 — ?■ L + 1/2 (that is, 5 — )> 0) as follows. 

Q.-i(x,y)--l^^^^_^ ((^ + i)/2)„.! ^^^l 1/2-. '1-^J 



/-L, 1/2 - L - e o^ 

^^\ 1/2-L '1-^J 



11 (1/2 -L -5)^(1/2), 
2 5 ((d+l)/2)^L! 

1 f^ 1 (l/2-L-e)Jl/2)„ /-n,l/2-L-e., A 

2^^_^^n-L-e {{d + l)/2)^nl ''\ 1/2 - n '^ ^J 
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Clearly (after shifting the index in the infinite series above) 

hm/o rxvl ll (l/2-L-5)Jl/2)^ (-LM2-L-e A 

2^L-n ((rf + l)/2)„n! ^ \ 1/2-n ' J 

1 - 1 (1/2 - L)„^^(l/2),^^ /-L-n,l/2-L , 2\ 

2^n((rf+l)/2)„^,(n + L)!^ ^^ 1/2-L-n ' ^/ 

We remark that the hypergeometric polynomial 

/n(.):=2Fi^ 1/2-L-n ' ^ " ^' J 

is strictly monotonically increasing on (0, 1) and assumes there its maximum value 1 at 
V = 1. For example, this can be seen from applying linear transformations of variable 
P, Eq.s 15.8.1 and 15.8.6 15.6.1] and change to integral representation [U Eq. 15.6.1] as 
follows 

}n[.V)-V 2-1^1 1 /9_ r_^' 1-7^ 



1/2 - L - n' y2 
' —L — n 



n ^ 



[1/2- L-n)^ 



2^V L + 1 ' 



-L-n)„ r(L + l) ^,, f' ty'-H^-tf^'''^\ ^ 



(1/2 - L - n)„ r(l/2) r(L + 1/2) io (1 - v' t) 

and the fact that the leading coefficients are positive. Hence, the series in the last limit is 
uniformly convergent with respect to v in [0, 1], since L > 1. 
Putting everything together, we arrive at 

/Ci+i/2(x,y)=lmi{Q^_i(x,y) + (-l)^+ic,(/3)||x-yf^-^} 

-limio fxvl ll (l/2-L-e),(l/2)^ f-L,l/2-L-e 
-hm|Q,_i(x,y)--- ((^ + i)/2),L! '^\ 1/2 - L '^ 

+ lim|iliiZ^-A^#i^.F/-^'3/2-f-;l_,A 
.^o|2e ((d + l)/2)^L! ' 'V 1/2 -iv ' y 

+ (-l)^+^c,(/3)||x-yf^-^ 



v^ 



ll(l/2-L-5)Jl/2), 



lim|Q,_,(x,y)--- ((^+i)/2)^^. 



2F.(-'J/1_-r^;i-.^)} 
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1(V2-L)^(V2)^ if f -L,s 1\ 

+ 2 ((rf+l)/2),L! ^ H-eY^\l/2-D^--2)- 

(54) + U. {^ f/;^;^-;^;;^^^- ^- + (-l)-c.(^) |,x - y I,-- 

In the last step we use a linear transformation of variable for the hypergeometric function. 
First, we observe that 



e^o^p ^\^l/2-L' ^;2y r e^o4^ (1/2 - L)^fc! ^ ^ 



tt ^ (V2-L), ^ ^ • 

By the definition of a Pochhammer symbol as a rising factorial, we have 

(1/2 -L-e)^ ^ /L + 1/2 + A ^ ^ (1/2 + e)^ ^ r(L + 1/2 + e) 

L! ^ ^ V ^ / L\ ^ ' r(l/2 + e)r(L + l) 

Furthermore, the duplication formula for the Gamma function gives 

.... /g , . _ 1 2^"^^ r((rf + l)/2) r(L + 1/2 + e) TjL + 1/2 + e) 

■^^ ^ ^ sin(7r£) 2 (L + £) r(L + (d + l)/2 + e) r(L + e) 

^ r(e)r(l-e) 2'-'" T{{d + 1)/2)T{L + 1/2 + e)T{L + 1/2 + e) 

71 2{L + e) T{L + {d + l)/2 + e)T{L + e) 

^ r(l + e) r(l - e) 2^-^^ T{{d + l)/2) T{L + 1/2 + e) T{L + 1/2 + e) 
£1(1/2)1(1/2) 2 r(L + (d + l)/2 + £)r(L + l + £) 

Hence, the expression in curly braces in ( 154|) can be rewritten as follows 



1 (-1)^ r((rf + l)/2) TjL + 1/2) r(L + 1/2 + e) ^^ 

e 2 r(l/2)r(L + (c/+l)/2)r(l/2 + £)r(L + l) 
1 (-1)^+1 r(l + e) r(l - e) T{{d + l)/2) r(L + 1/2 + e) T{L + 1/2 + e) ^^^,^ 
e 2 r(l/2)r(l/2) r(L + (t/ + l)/2 + £)r(L + l + £) 

(-1)^+1 T{{d + l)/2) T{L + 1/2) r(L + 1/2 + e) ^l ^(5) - 1 



2 r(l/2) r(L + (rf + l)/2) r(l/2 + e) T{L + 1) 
where the function 



r(i/2 + £) r(L + i) r(L + (d + i)/2) r(L + i/2 + £)r(i + £)r(i-£) 2, 
^^''~ r(i/2) r(L + i + £)r(L + (d + i)/2 + £) r(L + i/2) r(i) r(i) "" 
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assumes the value 1 at e = and is differentiable at with (using product rule) 

G{e) - 1 



G'(0) = lim 



T'{l/2 + e) 



r(l/2) 



r(L + i)r'(L + i + £) 



e^O 



[T{L + l + sW 



T{L + {d+ l)/2) T'{L + {d+ l)/2 + e) 



[T{L + {d + l)/2 + e)f 



+ 



T'{L + 1/2 + e) 



T'{l+e) 



e^O 



Til] 



e^O 



T'(l-e) 



e^O 



Til] 



e->0 



r(L + i/2) 

+ lnw^ 

Using the Digamma function ip{z) = r'{z)/ r{z), we can write 

G"(0) = ^(1/2) - ^(L + 1) - ^(L + {d+ l)/2) + ij{L + 1/2) + ^^(1) - ^(1) + \nv^ 
= 2 ln(2t;) - 2 In 2 + ^(1/2) + i){L + 1/2) - ip{L + 1) - ^(L + {d + l)/2). 

Thus, we obtain 



^.^.lljy2-L-eUimL^,, 



-l)^+iQ(/3)||x-y| 



£^o|2e ((rf+l)/2)^L! 

_ (-1)^+^ T{{d + l)/2) r(L + 1/2) r(L + 1/2) 

2 r(i/2)r(L + (d + i)/2)r(i/2)r(L + i) 

_ (-1)^+1 (l/2),(l/2)^ ^,, 



t;'^G"(0) 



(55) 



2 ((o? + l)/2)^L! 
-^(L + l)-V^(L + (t/+l)/2)). 

Putting everything together, we have 



w^^ 21n(2tO - 2 In 2 + ^(1/2) + ip{L + 1/2) 



r f \ r !n f \ 1 1 (1/2 -L- e) Jl/2)^ ^ /-L, 1/2 - L - g , 
/C,^,/,(x,y)=hm|Q,_,(x,y)---^^^^-^^^^^-^^,F,^ ^ /^ _ ^ ;1 



+ 



l (l/2-L)Jl/2)^ 
2 ((rf+l)/2)^L! 



-"Sli^^^Vl/2-L 



^'"■l-ll-l 



+ lim 



ll(l/2-L-e)Jl/2), 



e^o|25 ((d + l)/2)^L! 



1 ^"' 



1 (l/2-L)„(l/2)^ 
2^^-^ ((rf + l)/2)„d 



^^^ + (_l)^+ic,(/3)||x-y 



/-n, 1/2 - L. 2^ 

^^^l^ 1/2-n '1-^J 



1 2/3-1 



(56) 
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2{r{n{{d + l)/2)^^,{n + L)\' \ 1/2-L-n ' ) 

l {l/2-L)^{l/2), ' (-1)^ i-L), ,(,_,, . _ ,.k 
^2 {{d+l)/2),L\ ^ k (1/2 -L), ^ > 

+ ^~V ((rf+\')%^,L! ^'''(21^(2^0 - 2 In 2 + ^(1/2) + ^(L + 1/2) 
-^(L + l)-V^(L + (rf+l)/2)). 

The series expansion of Q_l_i/2 follows. 

Asymptotic analysis. Inspection of the derivation of the asymptotic relation in The- 
orem [H] shows that one is permitted to take the limit as s — )■ L + d/2 (which corresponds to 
2/3 — 1 — 7- 2L. In particular, the function Vt{k) remains uniformly bounded for sufficiently 
small e {s = L + d/2 + e) and therefore, by ( HHI) 



Af{L+d/2,d)^-2''-^ 



T{{d + l)/2) r(L + {d+ l)/2) r(2L + d) 
V{l/2-L)V{L + l + d/2) 



2 



T{d/2) 



-As 



V{d/2 + 2L + d) 



This completes the proof. D 

Proof of Theorem [II From the arguments following the definition of the kernel ([7]) we al- 
ready know that /C/3(x, y) is a symmetric positive definite kernel; that is the expansion 

oo 

(57) /C;3(x,y) = ^A„Z(rf,n)Pi'^)(x-y). 

n=0 

has positive coefficients Ao,Ai,... (cf. I. J. Schoenberg [20]). By Theorems [8] and [9] and 
the remarks following these theorems, the dominant term of A„ as ?i — )■ oo comes from the 
coefficient in the ultraspherical expansion of /C^(x, y) — Q^_i(x, y); that is, A„ x n~^^, 
where s = (3 + {d - l)/2. 

The Sobolev space H^(§'^) defined by the coefficients Aq, Ai, . . . , in (IF?]) is a reproducing 
kernel Hilbert space for s > d/2 (that is, /3 > 1/2) with reproducing kernel /C^(x, y). 
By uniqueness of the reproducing kernel Hilbert space 'H^(/C^,S'^) (uniquely defined by 
]Cp{yi, y)), n^iJC/s, S'^) and ^{8"^) coincide. D 

Proof of PropositionlE' Let d>2 and /3 > 0. By the Funk-Hecke formula 
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The standard substitution {1 — t)v = u — t gives 



(y -z 






(1-t) 






dv 

d/2-l 



a + ty 



■dv. 



-^ Jo {l + ]^,vY~'^' 

By [H Eq. 15.6.1] the last integral represents a Gauss hypergeometric function; that is, 

ff .N/3-1, .s ., (s+d/2~i d/2^iUd^iI\l3)T{d/2) fl-d/2,^ l-t\ 

JJ^-'-'^-^ da.(y) = (l-t) (1 + t)/ -^ ri^^d/2) '^\p + d/2-TT-t)- 

A linear transformation of variable (cf. [H Eq. 15.8.1]) yields 

(y . z - tt^ da.(y) = (1 - t)^-^/-^ ^^^^^^^2^/-^ .F, A " ^/2, ^/2. 1^ 



ood T{(3 + d/2)^ ''^\ (3 + d/2 
Using dl]) and simplification gives the result. D 

Proof of Proposition^^ Let d >2 and /3 > 1/2. Substituting the integral representation 
([7]) into the double integral and interchanging integrals, we get 

/ / /C^(x,y) dcrd(x)d(Td(y) 

= [ [ [ / (x-z-t)J-^(yz-t)^-Marf(z)dtdarf(x)darf(y) 
J§d >/§<* J-i Jst* 



= [ [ \[ (x-z-t)^-Mcr,(x)l [/ (yz-t)^-Ma,(y) 
j-i Js'i Us'i J L^'S'' 

By Proposition [2] the integrals in square-brackets do not depend on z. Hence 

/C/3(x,y) dord(x)d(Td(y) 

1 _ j.\W+d~2 



dcrd(z)dt. 



-,d-2 



r((rf + i)/2)r(/3)]' f\^ ,,2/3+d-2r ^ A-rf/2,rf/2. 1-t 



/?r(/3 + rf/2) _ 
The standard change of variable 2x = 1 — t gives 

/C^(x,y) dcrd(x)do-d(y) 



2F1 



/3 + d/2 



dt. 



-)d-2 



■r((rf + i)/2)r(/3) 



-\ 2 



22/3+d-i , ^ 
'0 



2/3+d-2 



A-rf/2,rfA; 



H 2 



/^ r(/3 + d/2) 

For d = 2 the hypergeometric function above reduces to 1 and therefore 



dx. 



/C/3(x,y) do-2(x)do-2(y) 



• r(i + i/2)r(/3) 
. Ar(/3 + i) 



Jo 



-\ 2 



2/3 



22/3+1 

2/3 + 1' 
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For d > 2 we apply Proposition [T7] to obtain 

1 2 



41 



sd J n't 



/C/3(x,y) d(Td(x)do-d(y) =2 



d-2 



r((rf + i)/2)r(/3) 

7rr(/3 + rf/2) 



g2g+d-i r(2/3 + d-i)r(i) 
r(r(2/3 + rf) 






2/3 + d - 1 
2/3 + d 



l-d/2,d/2 
P + d/2 



1 - d/2, d/2 
/3 + d/2 



1,1 



For d an even positive integer, the series expansion of the given Kampe de Feriet function 
terminates after finitely many terms (cf. ([22]) )• The result follows. D 



Proof of Proposition From fH9|) we have that 



QM-i(x,y) 



2Af - 1 ^-^ 



1 - M) Jl/2 - M)„ 



„=o ((rf+l)/2)„r2! 



-n. 



1-M 1-x-y 



2^11 3/2-M ' 2 



Expanding the hypergeometric polynomial and using that {—n)^/n\ = (— l)™/(ra — rra)!, 
we obtain after rearranging terms 



QM-i(x,y) 



1 "^^ (-1)"^( 
2M-1 ^ (3/2- 

m=0 ^ ' 



1-MX 



(3/2 -M)^m! 



'M-l 



(l-M)„(l/2-M)^ 
((d+l)/2)>-m) 



1-x-y 



The square-bracketed expression is a hypergeometric polynomial evaluated at unity and 
can be reduced as follows (cf. [H Eq. 15.4.20]): 



:^-MUl/2-M)^ ,m 

— 2-ri 



1 -M,m + l/2-M 
m + {d+l)/2 



(1 - M)„(l/2 - M)„ r(m + {d + l)/2) r(2M - m + d/2 - 1) 



i{d + l)/2\ 



T{M -l + {d+ l)/2) T{M + d/2) 



_^.m^2M-i (^/2)2M-i (1 - M)^(l/2 - M\ 



(d) 



2A/-1 



(2 - rf/2 - 2M), 



Thus 

(58) 

QM-i(x,y) 
(59) 



i2A/-l 



Af-l 



(^/2W-i >R (l-M)^(l-M)„(l/2-Af)^ A-x.y 



2M - 1 (d) 



E 



2Af-l ^=0 

2A/-1 T7 



(3/2-M)^(2-d/2-2M)^m! 



2^^''-^ (c?/2)2M_i p A - M, 1/2 - M, 1 - M 1-x-y 
2M-1 {d)^,,_, '''\3/2-M,2-d/2-2M^ 2 



The constant Cd{M) can be obtained from Theorem [81 



42 J. S. BRAUCHART AND J. DICK 

Since Qm-i(x, y) is a zonal function (depending on x-y), using (l59l) . by the Funk-Hecke 
formula and subsequent change of variable 2v = 1 — u we get 



/ QA/-i(x,y)d(Td(x) 






2M-1 (rf)^,,. 

^2M^ (rf/2)^ a;,_i /-^ ,/2_i ,/2-i p A - M, 1/2 - M, 1 - M. \ 

2M - 1 (rf),,,„, ' T;7yo ^^^ 3F,^3/2_M,2-rf/2-2M'^J^^' 

(^/2)2A.f-i od^i ^d-ir(rf/2)r(ci/2) A-M,l/2-M,l-M,rf/2 A 
1 (d)2,,„i cud T{d) ''\3/2-M,2- d/2 - 2M, rf ' ^ " 



22A/-1 

2M 



The last step follows from the terminating series expansion of the hypergeometric polyno- 
mial or [T^ Eq. 7.2.3.9]. Using ([T]) and simplification gives 

fn ( \A ( \ 2^^-^ (^/2)2M-i ^ [l-MM2-M,l-M,d/2 ,\ 
y^^QM-i(x,y)da.(x) = ^^^_^ (d),,,„, ^^3^3/2_,^^,2-^/2-2Mi'V- 

The 4 Fa-hypergeometric is Saalschutzian (that is, 1 + the sum of the upper parame- 
ters equals he sum of the lower parameters). Application of Whipple's transformation [H 
Eq. 16.4.14] yields 



A-Af,l/2-M,l-M,rf/2 \ (l)A/-i(3/2-rf/2-M) 

^^\d, 3/2 - Af, 2 - d/2 - 2Af '^) (2 - d/2 - 2M)^,^,{3/2 - 



X 4-^3 



P A - M, 1/2 - M, M + d - 1, d/2 \ 
^■'' d,l-M,(rf+l)/2 '^'• 



(Note that one must not simplify the formula by "cancelling" 1 — Af in the hypergeometric 
polynomial (evaluated at 1) at the right-hand side.) From 

22^^-1 {d/2),,,_, (l),,_,(3/2 - d/2 - M)^_, ^_^^,,_^ (A/ - 1)! 



2Af-l {d\,,^, (2-rf/2-2Af),,_,(3/2-Af)^_, ^ ' (3/2),,„i 

we obtain the result. 

For d = 2, we get with the help of Mathematica 

XJs.^""^^"''^^^"^^"^^"^^'^^^^"'^ W^)Z/\ 1-M,3/2,2 ' 



-1) 



^,„i (M - 1)! ""^ (1/2 - M)SM + 1)„ 



^E 



(3/2)m-i t^ (3/2)„(2)„ 

22M r(i/2- Ar)r(M) 



Af2(2Af + l)^ M(2Ar + l)^y? 

o2A/ ^_i)M-i (M-1)! 



+ 



Af2(2Af + l)^ Af(2Ar + l)^ (1/2) 



M 
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D 

Appendix A. The reproducing kernel Hilbert space ■H^(/C/3,§'^) 

In this section we show that all functions in ^{^(/C^, S^) have an integral representation 
of the form ([8]). Note that the set of functions Yl^=i '^i^pi.'^ Yi) is dense in "H^. Let / G T^^s. 
Then there exists a sequence of functions /i, /2, . . . G "H^ with 

n 

such that 11/ — fnWiCfi — )■ as n — )■ oo. Further let 

n 

Thus 

/„(x) = / / 5(„(z,t)(z-x-t)'^"Mcrd(z)dt. 

Then elementary algebra shows that 

Wfn- fm\\l, = {fn- fmJn- fm)Kp= / Ifi-nlz, t) - 5f„(z, t) |^ dad(z) dt. 

Since the sequence (/„) converges in the Hilbert space Tip, it is a Cauchy sequence. Thus 
ll/n — /m||iC;3 — >■ as n, m — )■ cxD and therefore the sequence (gn) is a Cauchy sequence in 
the Hilbert space L2(E>'^ x [—1, 1]). Thus the sequence (gn) converges to a function g in 
L2(§'^x [-1,1]). Let 

/(x)=/" /" (7(z,t)(z-x-t)^-Ma,(z)dt. 

Then we have 

o<ii/-7ik, <ll/-/„lk, + ll/„-7lk,. 

We have ||/ — /„|k^ — > as n — ?> oo by the definition of the fn and 

ll/n-7llL= / / k(z,t)-^(z,t)|'da,(z)dt^O 
J-lJ§d 

as ?T. — )■ cxD, since g is defined as the limit of the sequence (gn) in L2(S'^ x [—1, 1]). Thus we 
have 11/ — fWjCfj = 0, which shows the statement. 
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Appendix B. Auxiliary results 
Proposition 14. Let \a/A\ + \b/A\ < 1. Then 



"^2' 2/3,2/3" AM 



2A 



2A - a - 6 



a/2,(a + l)/2 

f3 + l/2,f3' + l/T\2A-a-b 



2A-a-b 



Proof. Since \a/A\ + |&/A| < 1, one can use the series expansion of the Appell function. 
Reordering gives 



m=0 n=0 

oo 



fia/A,b/A) = J2J2 

m=0 

oo 

= E 

m=0 

oo 



(")m+n(/^)m(/3'),„ fa\^ fb 



oo 



(2/3)„,(2/3')„m!n! VA^ \A 



m\' /m 



n=0 



A 



(2/3')„^! 



(2/3)„. m! 

(2/9L^!' 'V 2/3' MJU 






m=0 ^ '^^"* 

Since | ph(l — b/A)\ < tt, the quadratic transformation of variable [U Eq. 15.8.13] can be 
applied; that is 



f{a/A,b/A) 



2A 



y- («)m(/^)m ^p^ { {m + a)/2, {m + a + l)/2. 



2A-6y ^^ (2/3) m! 

m=0 ^ ' ''"T' 



(3' + 1/2 



2A-6 






Invoking the series expansion of the hypergeometric function and identities for the Pochham- 
mer symbol, one gets 

it I A I. / A\ oo oo ^ ^ ^ ^^ ^^ N /ON // , , 1 N /r>\ / o \ m / , \2?1 

/(a/A,&/A) - --- - --- - 



2A \" 



\2A-b) 



EE 

m=0 n=0 

oo oo 



i^U^Ui^ + «)/2)n.((rn + a+ l)/2), 
(2/3)„(/3' + l/2)„m!n! 



V-y-(«L(/3L2-2"("^ + «) 



m=0 

oo oo 

EE 



^ (2/3)„(/3' + l/2)„m!n! \2A-bJ \2A - b 



2a 



2a 
2A-b) \2A-b 

2n 



(")m+2n(/5) 



+2n\r^/rn 



m=0 n- 

oo 

n=0 



- (2/3)^(/3' + l/2)„ m!n! V2A - 6; \2A-b 



2a 



b/2 



[a. 



2n 



W' + l/2)^n\ 



2F1 



Application of [H Eq. 15.8.13] gives 



2n + a,/3 2a \ / 6/2 
2/3 '2A-bJ \2A-b 

(")2n 



2n 



2n 



2^-a-V ^(^' + l/2)„n! 

A2n + a)/2,(2n + a + l)/2. 
"" 2^M /3 + 1/2 



2a 



4A - 2a - 2a 



6/2 



2A - a - 6 



2n 
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Reversing to series expansions and invoking properties of the Pochhammer symbol, we 
arrive at 

f{a/Ab/A) _^^ («),^^,„ / a/2 V^ ( b/2 ^ ^" 

(2^^)" ;^o^(/^ + V2L(/3' + l/2Lm!n! {2A-a-bJ ^-a-b^ 

00 00 1^ ir,\ //„ I i\ /r,\ / \ 2m / 7 \ 2n 



yy («/2L+J(a + l)/2)„.,„ 



^=o^(/^ + V2L(/3' + l/2Lm!n! \2A-a-bJ \2A-a-b^ 
By [H Eq. 16.13.4] the double sum is an Appel F4-function. The result follows. D 

Next, we expand a hypergeometric polynomial in terms of normalized Gegenbauer poly- 

nomialsPf(t) = Cr)/^(t)/Cr-)/^(l). 

Proposition 15. Let c be not an integer < 0. Then 

where 

a, _ (-lf(-n),ib),{d/2) 



p fk-n,k + b,k + d/2 A 
^\ k + c,2k + d 'V' 



Z{d,k) {c)kid)2k 

Proof. The polynomial Qn{t) has degree n and can therefore be expressed in terms of the 
family of normalized Gegenbauer polynomials (ultraspherical polynomials) {Pq , . . . , Pn } 

n n 

Quit) = Y.^kPlt\t) = E zrZIT^^^'^)^^^'^^'^' 

fc=0 fc=0 ^ ' ' 

where the coefficients a^ depend on Q„ and 



l\n{t)Pi^'\t){l-t'f'^'dt. 



Z{d, k) LOd 

(Note that ^ /^^ Plf\t)P^^'^\t){l - t2)'^/2-i dt = 1/Z{d, k); cf., for example, p.) 
Using Rodrigues' formulas (cf. [H Eq. 18.5.5]) 

and /c-times successive integration by parts, we arrive at 



Z{d,k) Ud 2\dl2)J_^ 



J_^Qn{t){{l-t^r'^'^'fdt 



"^'-^ ^ l"Q^^\t){l-t')'^'^'-'dt. 



Ud 2>^{d/2)J_, 
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The differentiation rule for hypergeometric functions [1] Eq. 15.5.1] yields 
qWu) = (-l)M-^)fcWfc ^p^ fk-n,k + b. 1^ 



2^= {c)i^ ' \ k + c '2 

Using the (terminating) series expansion of Qn , we get 

Z{d,k) u, 2^{d/2\ 2^ (c), ^ {k + c)^v\ ""' 

where the standard substitution 2x = 1 — t leads to a beta integral 

= 2"^+"-^ f x'^+'^+'^/^-l (^ _ ^^k+d/2-l^^ 

Jo 

2,+d^ini^ + k + d/2)T{k + d/2) 



2 



r{u + 2k + d) 



^ ^j,^,_, T{k + d/2) Tid/2) {k + d/2)^{d/2), 

T{2k + d) {2k + d)^ 

^ ,_, V{d/2) r(rf/2) , {d/2), jk + d/2Ud/2), 
T{d) {d),, {2k + d)^ ■ 

Hence 

ak ^ ,j^i OOd-i T{d/2) T{d/2) {-l)\-n),{h),{d/2\ "^ {k - n)^{k + b)^{k + d/2)^ 
Z{d,k) u, r{d) {c),{d),, ^ ^k + c)^{2k + d),u\ 

_ {-in-n),{b),{d/2). 



{c),{d) 



kK"'J2k 



p fk-n,k + b,k + d/2 A 

^^\ k + c,2k + d ' y- 



D 



As a corollary to the last result we have the following 
Proposition 16. Let 3/2 — 13 be not an integer < 0. Then 



>F<-3/2'--/;^)-i:^^w*)^rw 



fc=0 

where 

^k.n 



Z{d,k) (3/2 



n\{l-l3),{d/2), fk-n,k + l-(3,k + d/2 \ 

2-P),{d),,{n-k)\' \ k + 3/2-/3,2k + d ' / 
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Letk + 1- (3>0. Then 

ak,n r(3/2 - /3) r(rf) n\ 



Z{d, k) r(l - (3) r(l/2) T{d/2) T{d/2) {d/2),{n - k)\ 



JO JO 

Proof. On observing that (—!)''"(— n)^ = n\/{n — k)\, the first form of the coefficients 
ak^n/Z{d,k) follows from Proposition [T5l 

Let k + 1 — f3 > 0. Then a hypergeometric function as in the coefficient ak^n/Z{d, k) can 
be expressed as follows (cf. [H Eq. 16.5.2] and p!| Eq. 15.6.1]). For b > and c > 

-n,b,c \ _ r(6 + l/2) f\b-i(^ /^l/2-l ■p/-n,c. 
v6 +1/2, 2c' Ij- 1(6)1(1/2)70 ^ ^ ^ ' '^ 2c ' 

r(6 + i/2) r(2c) f' f't'-'x'-'{i-tf^-'{i -xy-\ 

U. tXv U. (/ 



r(6)r(i/2)r(c)r(c)7o Jo (i-xty^ 

That is (with b = k + l — P,c=k + d/2, and n = n — k) 

r(A; + 3/2-/3) r(2A; + d) 



p fk-n,k + l- f3,k + d/2 \ 
^ \ k + 3/2-/3,2k + d 'J 



r(A; + 1 - 13) r(l/2) T{k + d/2) r(A; + d/2) 

Jo Jo 
The coefficient in front of the double integral equals 

r(3/2 - /3) Tjd) (3/2 - /3), (d),, 

r(i-/3)r(i/2)r(d/2)r(rf/2) (i-/3), {d/2),{d/2),- 

The second form of the coefficient ak,n/Z{d, k) follows. D 

Proposition 17. Let Re^y > Re /3 > and Re(c - a - b) > and Re(c' - a' - 6') > 0. 
Then 



a,b a\b' 



1,1 



/' t^- (1 - 1)-^- .F. (V; t) .F, («'/; t) d t = r(^)r(7-/3) ^g. 

Proof. Under the given assumptions the hypergeometric functions have series expansions 
that are uniformly and absolutely convergent for all x with |x| < 1. Let A denote the 
integral we want to compute. Substituting these series expansions and interchanging sum- 
mation with integration, we get 

m=On=0 y^)m\C)n^-^- JO 

Evaluating the beta integral and using Pochhammer symbols, we obtain 

^^ _ r(/3) r(7 - /3) y. y. m^^n JaUb)^ ia'Ub\ l^l" 

r(7) £^o^(^Wn (^)- (^')n '^■''■' 
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The double series above represents a Kampe de Feriet function evaluated at unity argu- 
ments. The result follows using f E2]) . D 

Appendix C. Limits 

Because of fIFIl) limit and integration can be interchanged in the first integral in (152 p . 
For the second integral in f l52|) we observe that the function G{x] e) has the following 
Taylor expansion with integral remainder 

oe e=o Jo (y^ ^'^^^ 

where _ 

— g^ — = (logx) +{logx)gi{s)+g2{e) 

and the functions gi{e) and g2{£) are continuous on [0, 1] and depend only on L. This 
shows that in the second integral in ( l52l) logarithmic singularities appear; that is, apart 

from the term with — q^^\ _q (that contributes to the desired limit) the remainder is 
1 y _2i-2*(l-/;)^ (x.z-yzf'- / / x-z-yz "" 



2A. 2 V. (2-x.z-y.zr«-^V'°H2-x.z-y.zl ' "^■'^'> + 

Since 1 — x-z — y-z>c>0fory7^x and y ^ — x and 2/3 > 2L + 1, the integrals above 
can be bounded independently of e. Thus, the remainder above vanishes as e — ?■ 0. Hence, 
we can interchange the limit in the second integral in (l52l) . 
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